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PREFACE 

Tim  short  course  on  laser  effects  was  developed  for  a  12-hour 
series  of  lectures  delivered  by  the  author  at  the  Naval  Postgraduate 
School  in  the  fa'I  of  1973.  The  lectures  were  part  of  Professor  John 
Neighbours’  course  on  solid-state  physics,  which  stressed  topics  ap¬ 
propriate  to  an  understanding  of  laser  effects  by  solid  state  physicists. 

The  author  is  grateful  to  Professo--  Neighbours  for  the  oppor- 

i-Ullliy  lu  li.Oatr  irCLuirs  <r.tiu  fui  iiuiiiy  Vcuiim/ic  Li iSCLLSolOiih 

during  their  development.  He  also  wishes  to  thank  Professor  Otto 
Heinz,  chairman,  and  those  other  members  of  the  Department  of 
Physics  and  Chemistry  of  tin  Naval  Postgraduate  School,  too 
numerous  to  mention  here,  who  through  their  support  and  hospi¬ 
tality  made  valuable  contributions  to  this  work. 


RESPONSE  OF  MATERIALS  TO  LASER  RADIATION: 
A  SHORT  COURSE 


1.  LASERS 
1.1.  Introduction 

The  word  laser,  of  course,  is  an  acronym  for  “light  amplification  by  the  stimulated 
emission  of  radiation,”  but  that  is  not  terribly  enlightening.  More  correctly  described  a 
laser  is  a  device  for  producing  light  that  is  almost  totally  coherent.  It  works  in  principle 
like  this:  An  atom  emits  a  photon  of  light  when  it  decays  from  an  excited  energy  state 
to  a  lower  state;  the  difference  in  energy  between  the  two  states  AE  determines  frequency 
v  according  to 


AE  =  hu  (1) 

where  h  is  Planck’s  constant.  This  is  illustrated  in  Fig.  1. 
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Fig.  1  —Energy  levels 


This  is  the  case  for  any  light  source,  whether  laser,  flame,  incandescent  body,  etc.  In 
the  conventional  light  source,  atoms  emit  photons  in  a  random,  sporadic  manner  and 
spontaneously  decay  to  lower  states  when  excited  by  heat  or  electric  current.  In  a  laser, 
on  the  other  hand,  the  photons  are  emitted  in  phase  and  the  electromagnetic  radiation 
thus  produced  is,  more  or  less,  simply  a  propagating  sinusoidal  radiation  field  that  can  be 
described  on  a  macroscopic  level  by,  for  example, 


&  = 


Re 


e-2nkz/\eiu)(t-nz/c) 


(2) 


where 


Note:  Manuscript  submitted  January  31,  1974. 
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t’  is  the  electric  field  of  the  radiation 

Re  stands  for  the  real  part  of  the  complex  quantity  in  brackets 
£>0  is  the  maximum  amplitude 
h  is  the  extinction  coefficient;  in  a  vacuum,  k  =  0 
z  is  the  direction  in  which  the  wave  is  propagating 
X  is  the  wavelength 
t  is  time 

n  is  the  index  of  refraction;  in  a  vacuum,  n  =  1 
c  is  the  velocity  of  light  in  vacuum. 

Equation  (2)  is  a  standard  representation  of  the  electric  field  of  a  traveling  light 
wave.  However,  if  one  measures  the  electric  field  at  some  point  in  space  for  light  from  a 
conventional  source,  the  sinusoidal  variation  expressed  in  Eq.  (2)  does  not  appear,  for  the 
atoms  emitting  the  light  are  doing  so  at  rrnd  '  ft*,  and  the  sinusoidal  variation  due  to  the 
emission  from  each  atom  is  averaged  to  some  .  ime-independent  value.  This  is  not  true  of 
laser  emission,  where  the  individual  photons  art  in  phase.  Measuring  the  electric  field  at 
a  point  in  space  for  laser  light  results  in  the  oscillating  £  predicted  by  Eq.  (2). 

This  coherence  is  created  by  taking  advantage  of  stimulated  emission  in  materials  in 
which  metastable  states  can  be  induced.  The  lifetime  of  an  atom  in  an  excited  state  de¬ 
pends  on  the  quantum  mechanical  selection  rules  for  transition  to  a  lower  state,  and  there 
are  states  from  which  transition  to  a  lower  level  is  extremely  improbable.  Such  states  are 
called  metastable,  and  an  atom  not  disturbed  by  outside  influence  will  remain  in  a  meta¬ 
stable  state  for  a  very  long  time.  If  a  metastable  atom  interacts  with  a  photon  of  fre¬ 
quency  such  that  A E  =  hv,  where  A E  is  the  energy  difference  between  the  atom’s  normal 
and  metastable  states,  stimulated  emission  will  occur.  The  atom  will  decay  to  its  normal 
state  by  emitting  another  photon  of  frequency  v,  so  that  the  net  result  is  two  photons, 
and  the  second  photon  will  have  the  same  phase  temporally  and  spatially  as  the  first. 

In  a  laser,  then,  one  establishes  a  large  number  of  atoms  in  metastable  states  and 
arranges  the  optics  to  increase  the  likelihood  of  stimulated  emission.  Schematically,  a 
typical  laser  oscillator  looks  like  Fig.  2.  The  pumping  radiation  (for  example,  light  from 
a  flashlamp)  excites  the  atoms  in  the  lasing  medium  (for  example,  Cr+  +  +  ions  in  ruby). 

In  the  decay  process  (if  we  have  a  successful  laser),  a  large  number  of  ions  are  left  in  a 
metastable  state;  this  is  called  a  population  inversion.  As  some  atoms  begin  to  decay, 
they  stimulate  others  to  decay.  But  this  alone  would  not  provide  a  laser,  since  the  emis¬ 
sion  would  occur  in  random  directions.  The  role  of  reflection  is  very  important;  the 
photons  moving  perpendicular  to  the  reflectors  pass  through  the  medium  many  times  and 
on  each  pass  more  and  more  atoms  are  caused  to  emit.  This  results  in  the  build-up  of  a 
very  strong  coherent  light  signal  that  travels  in  a  single  direction.  Useful  light  output  is 
obtained  by  making  one  of  the  mirrors  a  partial  reflector. 
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Fig.  2— Schematic  representation  of  a  laser 


It  is  interesting  to  look  at  a  few  examples  of  the  intensity  of  laser  light.  In  a  typi¬ 
cal  ruby  laser,  the  concentration  [1]  of  Cr't'  +  +  ions  is  about  2X  1019  cm-3,  and  popula¬ 
tion  inversions  are  of  the  order  of  3  X  1016  cm-3.  Crudely  speaking,  we  can  think  of 
creating  3  X  1016  quanta/cm3  in  the  lasing  medium.  Since  we  have  arranged  the  laser  so 
that  the  output  is  in  a  single  direction,  and  since  photons  move  with  the  speed  of  light, 
we  obtain  3  X  1016  X  3  X  1010  =  9  X  1026  quanta/cm2s  from  the  laser.  For  ruby,  the 
lasing  wavelength  is  6943A,  and  since  the  energy  of  each  quanta  is  hv ,  one  can  readily 
calculate  that  the  output  is  about  2.5  X  108  VV/cm2. 

Let  us  compare  this  to  the  power  that  a  hot  body,  say  the  sun,  emits  at  the  same 
wavelength  with  a  similar  bandwidth.  This  can  be  calculated  by  the  use  of  Planck’s 
radiation  law, 


ij  -  1 

w  ~  7T2c3  e'»«/*T  ' 


(3) 


Uw  is  the  energy,  per  unit  volume  and  per  unit  bandwidth,  radiated  by  a  blackbody  at 
temperature  T ;  k  is  Boltzmann’s  constant.  The  radiation  leaves  the  black-body  source  at 
rate  c,  so  the  power  radiated  per  unit  area  of  the  source,  per  unit  bandwidth,  is 


_  cUw  _  nw3  1/4 
W  4  "  j)-2c2  ehu>!hT  _  1  ' 

If  we  use  the  sun’s  temperature  of  6000  K,  and  A  =  694 3A, 

Iw  *=  2X10*5  erg/cm2. 

For  the  ruby  laser,  a  typical  line  width  is  3A,  so  Aw  %  1.2 X  1012  s-1.  Thus  the 
power  density  at  the  source  is 

/  ~  2.5  X  107  erg/cm2  s  ~  2.5  W/cm2  . 

Thus  the  power  density  for  comparable  narrow-bandwidth,  nearly  single-frequency 
light  is  much  greater  at  a  laser  source  than  at  a  conventional  hot-body  source,  because 
laser  light  is  coherent. 
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1.2.  Propagation 

i  hi'  propagation  of  laser  light  through  the  atmosphere  poses  a  complex  problem  and 
it  will  not  be  discussed  here.  Suffice  it  to  say  I  hat,  as  anyone  who  has  driven  on  a  foggy 
night  certainly  realizes,  1.,.  l  is  certainly  scattered  in  the  atmosphere.  Lasers  of  high  power 
density  post-  even  more  difficult  propagation  problems  because  the  high  intensity  warms 
the  ;ur  and  creates  a  density  change  across  the  beam.  Tins  variation  in  density  refracts 
the  light  and  causes  beam  spreading,  or  “thermal  blooming.” 

Consider  briefly  the  propagation  or  laser  light  in  free  space  or  in  vacuum.  Under 
these  ideal  conditions,  the  only  change  in  the  power  density  is  due  to  simple  beam  diver¬ 
gence.  Since  the  typical  laser  emits  light  that  is  nearly  unidirectional,  the  beam  diver¬ 
gence  is  small.  In  fact,  one  feature  of  a  laser  is  that  the  divergence  is  nearly  at  the  dif¬ 
fraction  limit,  which  is  of  fhe  order  of  A/a,  where  a  is  the  diameter  of  the  output  aperture 
of  the  laser.  For  the  ruby  laser  discussed  above,  this  gives  a  divergence  angle  of 


0  * 


6913 x  10'* 

r 


7  X  10' 12  mracl 


for,  say,  a  1-cm  aperture.  In  practice,  one  needs  to  go  to  much  trouble  1o  realize  this 
limit  of  divergence,  but  it  has  been  done.  More  commonly,  an  “off-the-shelf”  ruby  laser 
might  have  a  beam  divergence  of  a  few  mrad. 


The  newcomer  to  lasers  has  usually  heard  about  diffraction-limited  beams  and  the 
consequent  extreme  directionality  of  laser  right  He  is  usually  surprised  to  discover  that 
at  long  distances  from  the  source  these  beams  have  power  densities  that  vary  as  the  recip¬ 
rocal  of  the  squan-  of  the  distance,  !iho  all  radiating  sources.  To  see  this,  consider  a 
source-  of  power  P  W,  diameter  a,  and  divergence  angle  0,  as  shown  in  Fig.  3. 


At  distance  r  from  the  source,  the  power  density  is 


J 


_ P _ 

--  (u  +  2r  tan  <1  )  ' 


or,  since  U  is  very  small  and  tan  0  ^  (I, 
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j  (a  +  2 rO)2 


of 

From  this  expression  it  is  apparent  that  for  large  distances,  such  that  2 r/a  »  1, 


7 ro2  4 r202 


nr202 

or,  since  0  =  X/a, 

r  =  _P_a£ 
nr2  X2  ' 

For  example,  consider  a  10-kW  beam  of  10.6-jum  wavelength  and  10-cm  aperture  at  1  mi 
(i.e.,  a  high-power  CO2  laser); 

7  =  _ 104  X  102 _ 

tt(5280X  12X  2.54)2(10X  10'4)2 


7=12  W/cm2. 

From  Eq.  (5),  if  we  substitute  Iq,  the  power  density  at  the  source,  for  7>/(7ra2/4)  and 
recall  that  0  ^  X/a, 


I  =  Io 


1  +  2r 


From  this  expression  we  can  see  that  if  r  is  small  there  is  little  change  in  power  density 
emitted  by  the  source.  The  distances  at  which  thir  s  true  are  referred  to  as  “near  field,” 
and  the  fine  details  of  the  beam  pattern,  such  as  local  variation  in  intensity,  hot  spots, 
etc.,  are  preserved  in  the  near  field.  It  is  apparent  from  Eq.  (6)  that  this  near-field  dis¬ 
tance  will  be  limited  to  r  such  that  /  =»  Iq,  or 


«  1 
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or 


rnear  field  A  • 

For  lasers  with  exceptionally  good  optics  that  have  a  Gaussian  distribution  of  power  den¬ 
sity  across  the  beam,  the  near  field  pattern  will  persist,  for  distances  on  the  order  of  a2j\ 
[2-4]. 

As  a  final  comment  on  power  densities  at  distances  from  laser  sources,  let  us  use 
Eq.  (6)  to  calculate  the  distance  at  which  the  power  density  is  halved: 

1  =  L 
A)  2 


and 

r,B  <>/r-D. 

For  our  illustration  of  a  CO2  laser  with  a  10-cm  aperture,  r  =  680  ft,  or  a  little  more 
than  0.1  mi. 


1.3.  Response 

The  remarks  here  were  not  intended  to  explore  lasers  or  laser  propagation  at  any¬ 
thing  beyond  the  most  basic  level.  They  were  designed  rather  to  set  the  stage  for  the 
main  purpose  of  these  notes,  that  is,  to  describe  the  behavior  of  ma tenuis  under  irradia¬ 
tion  by  laser  light.  A  great  deal  of  study  is  under  way  in  this  field,  for  the  application 
of  highly  intense  light  to  materials  has  many  uses.  In  the  medical  world,  for  example, 
lasers  are  used  to  “spot-weld"  detached  retinas.  The  Ford  Motor  Company  is  installing  a 
large  laser-computer  machine  to  automatically  weld  parts  of  Torino  bodies.  There  are,  of 
course,  possible  military  applications. 

We  address  ourselves,  then,  to  a  detailed  examination  of  what  happens  when  the 
light  interacts  with  the  material.  It  is  convenient  to  think  of  this  interaction  in  three 
parts,  although  all  three  influence  each  other  and  take  place  simultaneously. 

First,  the  light  couples  with  the  material.  To  a  first  approximation,  this  is  governed 
by  simple  optical  reflectivity,  although  at  high  power  densities  other  effects  become  im¬ 
portant.  The  net  result  of  the  coupling  is  the  conversion  of  some  fraction  of  the  optical 
energy  into  thermal  and/or  mechanical  energy. 

Second,  the  thermomechanical  signal  is  propagated  into  the  material.  The  details  of 
this  propagation  play  a  prominent  role  in  determining  the  net  effect  of  the  irradiation. 

For  example,  pure  copper  or  aluminum,  with  their  high  thermal  diffusivity,  can  readily 
dissipate  large  amounts  of  thermal  energy  and  thus  require  higher  intensities  for,  say, 
laser  welding  than  more  poorly  conducting  metals  like  stainless  steel. 
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Third  is  the  induced  effect,  the  effect  the  thermomechanical  energy  has  on  the  ma¬ 
terial,  such  as  melting  or  vaporization,  shock  loading,  crack  propagation,  and  so  on.  Here 
again  there  is  an  interplay  with  the  other  parts.  When  a  metal,  for  example,  begins  to 
melt,  its  optical  reflectivity  and  thermal  diffusivity  change  markedly,  and  this  changes  the 
coupling  and  the  energy  flow.  So  in  our  discussion  we  shall  have  to  be  aware  of  the  in¬ 
teraction  between  the  three  aspects  of  coupling,  energy  flow,  and  induced  effects. 


2.  OPTICAL  REFLECTIVITY 
2.1.  General  Properties 

To  consider  the  coupling  of  the  laser  energy  to  a  material,  we  need  first  to  know  the 
optical  reflectivity  R  and  the  transmissivity  T  for  light  incident  on  a  surface  which  divides 
two  semi-infinite  media.  The  transmissivity  plus  the  reflectivity  equals  unity  at  a  single 
surface : 


R  +  T  =  1 .  (7) 

In  most  practical  situations  we  are  dealing  with  more  Uv.n  one  surface;  typically,  we  have 
a  slab  of  material  with  light  impinging  on  one  surface.  Some  light  is  reflected,  and  the 
rest  is  either  absorbed  or  passes  completely  through  the  slab.  In  such  a  situation  we  shall 
describe  the  net  result  of  all  the  reflection,  after  multiple  passes  inside  the  slab  and  ap¬ 
propriate  absorption  has  been  accounted  for,  in  terms  of  the  reflectance  :R,  the  absorp- 
tance  U,  and  the  transmittance  f!  : 


-R  +  ,T  +  (?  =  1 .  '  (8) 

What  we  really  are  interested  in  from  the  point  of  view  of  material  response  is  (f,  the 
absorptance  of  the  material.  In  most  materials  of  interest  from  the  practical  aim  of  using 
lasers  to  melt,  weld,  etc.,  c)  is  zero,  and 

')?  +  (■?=  1.  (9) 

In  a  later  section  we  shall  consider  the  relationship  between  .'K  and  R. 

To  understand  reflectivity,  we  must  use  some  general  results  from  the  theory  of 
electromagnetic  waves.  Let  us  summarize  these  briefly  at  this  point.  The  electric  field 
of  the  electromagnetic  wave,  from  Eq.  (2),  is 

£  =  R  e[£,oe-2nkz/\eiu(t-nz!c)] 


The  relationships  we  need  are  those  among  the  index  of  refraction  n,  the  extinction  coef¬ 
ficient  /j,  and  the  material  properties.  These  relationships  can  be  derived  by  substituting 
Eq.  (2)  in  the  wave  equation 


<)z  2 


Me 


+  no 


(10a) 


This  results  in  the  expression 
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pc(- CO2)  +  tcopa. 


(10b) 


Note  that  we  are  using  rationalized  MKS  units  throughout.  The  material  properties  enter 
through  p,  e,  and  a,  which  are  the  magnetic  permeability,  the  dielectric  function,  and  the 
electric  conductivity  of  the  medium.  Using  the  usual  equations  between  the  field  vectors 


we  have 


KeCQ$  , 

(Ha) 

KmV  oH, 

(lib) 

a& , 

(11c) 

Kec o  . 

(lid) 

^771^0  • 

(lie) 

In  Eqs.  (11),  H o  are  the  electric  permittivity  and  magnetic  permeability  of  a 

vacuum.  Ke  is  the  dielectric  constant  and  Km  the  magnetic  permeability  of  the  material. 
By  substituting  Eqs.  (11a)  and  (lie)  into  Eq.  (10b)  and  using  2n/\  =  to/c,  we  obtain 

2 

(fe  +  w*)2  =  -KcKme0H0C2  +  iKmHoa  —• 

Finally,  if  we  introduce  c2  -  (eoPo)"1  anc*  ^o  some  algebra, 

n  -  ik  =  y/KZ  /Kc  -  i  J~-  (12) 


This  equation  relates  the  material  parameters  Km,  Ke,  and  a,  which  in  general  may  be 
complex,  to  index  of  refraction  n  and  extinction  coefficient  k.  To  describe  the  propaga¬ 
tion  of  the  light  wave  thus  requires  a  knowledge  of  Ke ,  Km>  and  o.  Before  we  describe 
these,  let  us  look  at  two  more  general  properties  of  our  propagating  electromagnetic 
wave. 


The  first  of  these  is  absorption.  If  the  medium  is  absorbing,  the  intensity  will  fall 
off  to  1/e  of  its  initial  value  in  a  distance  5,  obtained  by  setting  £2  of  Eq.  (2)  equal  to 

(l/e)frmaX.  or 


4nk5 

X 


1 


=  X 
4tt/j  ’ 


(13) 


This  shows  why  It  is  called  the  extinction  coefficient,  for  it  determines  skin  depth  X. 
Equation  (13)  is  fairly  general  in  that  once  k  is  known,  5  can  be  calculated.  As  noted, 
a  knowledge  of  the  material  properties  is  required  to  calculate  k. 
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The  second  general  property  we  wish  to  derive  is  the  expression  for  reflectivity,  in 
terms  of  n  and  k.  To  do  this,  consider  light  impinging  normally  onto  an  ideal  solid  sur¬ 
face,  as  shown  in  Fig.  4.  Here  we  have  illustrated  the  incident  (6,-),  reflected  (&r),  and 
transmitted  (£»*)  electric  waves  at  a  vacuum— material  interface.  For  the  present,  we  limit 
our  discussion  to  the  case  of  normal  incidence.  We  now  consider  the  boundary  condition. 
We  have 


&i  +  &r  =  &t  (14) 

for  the  electric  field.  For  the  magnetic  field  B,  we  write 

Bj  —  Br  —  B  t  .  (15) 


Fig.  4— Incident,  transmitted,  and  reflected  electric 
vectors  at  an  interface 


The  minus  sign  is  before  Br  because  &  X  B  is  positive  in  the  direction  of  propagation  of 
the  wave.  Now,  the  relationship  between  B  and  or,  since  B  =  juH,  between  H  and  &, 
is  required  in  order  to  proceed  further.  This  follows  directly  from  Maxwell’s  equations: 

V  X  £  =  -n  M  (16) 

VXH=  „6  +  e-f-  (17) 

It  is  convenient  to  rewrite  Eq.  (2)  and  introduce  cuX  =  27rc,  to  have  &  explicitly  in  terms 
of  a?  instead  of  both  u>  and  X.  Recall  that  &  is  a  vector,  and  take  it  as  being  along  the 
x  direction.  Thus 


/  co 


6*  = 


e 


z(n-ik) 


(18) 
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Here  we  have  dropped  the  “Re”  notation,  and  shall  simply  note  that  we  always  mean  the 
real  part  when  we  write  the  wave  in  exponential  form.  We  shall  use  unit  vectors  x,  y, 
and  z. 

Now  the  curl  expressions  reduce  to 


y  X  6  =  y 


bz 


which,  with  Eq.  (16),  tells  us  that  H  has  only  a  y  component,  so 

9HV 


7  X  H  =  -x 


bz 


Thus  Eqs.  (16)  and  (17)  become 


3  fc'x  bHy 


bz 


-M 


bt 


and 


bHy  „  X.  0t'x 

-7—  =  at'x  +  c  — r~ 
dz  bt 


(19) 


(20) 


and.  of  course,  &y  =  =  Hx  =  Hz  =  0.  Putting  the  expression  for  from  Eq.  (18) 

into  Eq.  (19)  leads  to 


I  CO 


Hy  = 


n  -  ik  c  --z(n-ih) 


HC 


f’liC  c 


This  is  the  desired  relationship: 


Hy  - 


(21) 


At  this  point  we  note  in  passing  that  Eq.  (20)  or  (10b)  could  be  used  to  yield  the  rela¬ 
tionship  of  n  and  k  to  m,  e,  and  o.  If  the  reader  is  unfamiliar  with  these  relationships,  it 
is  instructive  to  carry  out  the  algebra. 

Returning  to  our  consideration  of  the  reflected  electric  and  magnetic  fields,  we  re¬ 
write  Eqs.  (14)  and  (15)  with  the  help  of  the  relationship  between  H  and  <S,  from  Eq. 
(21); 


Sr  = 


C 


and 


l#l  -  A*l Hr  =  M2 Ht 


becomes 
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& . 


'll  ~  ‘M  s, 

”2  ~  /*i/ 


Solve  for  £,./£,•  by  eliminating 


£r  _  ni  ~  n2  ~  i(k i  ~  ^2) 
6/  ~  Hj  +  «2  -  '(fel  +*2) 


Finally,  the  reflectivity  R  at  the  surface  is 


R 


(«1  ~  H2)2  +  (fej  ~  k2)2 
(n  1  +  n2)2  +  (^i  +  ^2  )2 


(22) 


Take  medium  1  as  a  vacuum  and  drop  the  subscript  2.  This  gives,  since  in  a  vacuum 
ni  =  1  and  k\  =  0, 


(n  -  l)2  +  k2 
(n  +  l)2  +  k 2 


(23) 


Equation  (23)  is  the  second  relationship  we  will  ffhd  useful  in  discussing  the  coupling  of 
optical  radiation  with  metals.  Note  that  it  is  derived  for  the  special  case  of  normal  inci¬ 
dence  and  is  applicable  to  a  vacuum-material  interface. 


2.2.  Reflectivity  of  Metals  at  Infrared  Wavelengths 

We  turn  now  to  a  derivation  of  the  optical  reflectivity  of  metals  for  infrared  wave¬ 
lengths,  where  experiment  has  shown  that  the  free-electron  theory  (sometimes  called  the 
Drude-Lorentz  theory)  of  metals  is  adequate.  This  theory  rests  on  three  assumptions. 

The  first  is  that  electromagnetic  radiation  interacts  only  with  the  free  electrons  in  the 
metal.  The  second  is  that  the  free  electrons  obey  Ohm’s  law,  or,  more  specifically,  that 

*  dv  ^  m*  c-  /0^ 

m  dt  +  ~r~  u  ~  ~e ®  (24) 

where  m*  is  the  effective  mass  of  the  electron,  v  the  velocity,  r  the  relaxation  time  due 
to  collisions,  and  ~e£>  the  force  on  the  electron  due  to  the  electromagnetic  field.  The 
third  assumption  is  that  the  free  electrons  of  a  metal  can  be  described  in  terms  of  a  single 
effective  mass,  carrier  concentration,  and  relaxation  time.  There  has  been  a  good  deal  of 
discussion  about  the  validity  of  these  assumptions  in  the  literature.  Recent  work  [3]  in¬ 
dicates  that,  for  wavelengths  in  the  intermediate  infrared  (a  few  microns  to  many  tens  of 
microns)  and  beyond,  the  free-electron  theory  does  a  reasonable  job  of  predicting  the 
reflectivity  of  metals. 

To  derive  the  free-electron  optical  reflectivity,  we  try  solutions  to  Eq.  (24)  of  the 

form 


v  ~  eic°l , 


U 
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so  that 


and 


\  m* 

'"1  (<<-')  ^  ~ 


-  nt-, 


u 


i'T _ 

m*{  1  +  i  cor) 


6  . 


Now  the  current  flow  obeys 


o  -  ufe,  —  ~.\’cu 

whole  /V  is  the  electron  concentration  (number  of  electrons  per  unit  volume;.  By  com¬ 
parison  of  the  last  two  equations. 


o  er 

Ac  m  *  ( 1  +/cOT) 


or 


o 


_ iXe-T 

m  1 1  +  ituT) 


Now  the  dc  conductivity  is 


oo 


A !c2t 
rr,  * 


(25a) 


We  see  o  is  a  complex  quantity  and  seek  to  write  it  a;  the  sum  of  a  real  and  imaginary 
part.  Thus, 


iVe2r(l  -  icot) 
in*(  1  +  oi-r2  ) 

Define 

0  =  0,  -  1(12. 

The  result  is 

°u 

°]  ~  1  +  OJ2T2 

(25b) 


Oo 


(lowr 

1  +  t o2r2 


(25c) 


To  proceed  further  we  need  to  use  the  general  expression  for  electromagnetic  waves  de¬ 
veloped  in  Sec.  2.1.  Recall  Kq.  (12): 
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n 


v'AVi 


~~n~ 

CyOJ 


anu  from  tin*  complex  o, 


t/t 


Ol 

i  — 


-  to 
<-t)~ 


If  we  assume  only  free-electron  optical  interactions,  the  rnetal  does  not  polarize  under  the 
wave,  and  Kc  -1.  In  addition,  for  metals  in  the  infrared,  Km  =  1.  Thus, 


or 


n  -  ik 


Itj)  +  U  2 

fow 


n  ~  ik 


02 

r0^' 


«1 

loco 


(26) 


It  remains  only  to  separate  the  real  and  imaginary  parts  ol  Kq.  (26),  which  will  yield  two 
equations  in  n  and  It  and  thus  g>vo  n  and  It  in  terms  of  the  dc  conductivity  oL>  and  the 
rela'ation  time  r.  Then  we  can  use  our  expression  for  tiie  reflectivity  from  Eq.  (23)  to 
generate  H  from  it  and  h. 

To  carry  out  the  algebra  we  use  the  identity 


s/A  t  iB 


It  +  ,1 


!li  - 

V  2 


where 

It  -  y/X'1  +  It*  . 


Letting 


A 


1 


02 

tyco 


and 


II 


Q| 

C  0  OJ  ’ 


we  liave 


f 


(, .  i!L.f 

+  /£»  f 

\  <•« 

Vo^-V 

(27a) 


and 
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r  =  *)2  +  fe2 

(n  +  l)2  +  k2 

(Eq.  (23)),  give  the  reflectivity.  These  are  elaborate  expressions  indeed.  It  is  useful  to 
look  at  some  limiting  cases  and  at  the  material  parameters  that  determine  R. 

Notice  that  R  is  a  function  solely  of  Oj ,  o 2,  and  co.  Look  again  at  Eqs.  (25)  and 
note  that  Oq  can  be  used  to  replace  r  in  the  expression  for  a 1  and  o?: 


1  +  co2 


*2  2 
m**o  0 

N2c4 


(28a) 


a2 


1  ‘f  tJ2 


*2  2 
m*^o0 

N2e-1 


(28b) 


Equations  (28)  show  that  aj  and  02,  and  thus  /?, ‘depend  on  frequency  co,  constant 
m*/N ,  and  dc  conductivity  oq.  Thus 

R  -  f(oo.  a0,  m*/N) .  (29) 


This  means  that  we  can  use  the  dc  conductivity  to  predict  the  reflectivity.  Furthermore, 
if  we  know  the  temperature  variation  of  o 0,  we  can  use  this  method  to  calculate  R  as  a 
function  of  temperature.  This  is  a  useful  result,  because  it  is  difficult  to  measure  optical 
reflectivity  as  a  function  of  temperature,  whereas  it  is  fairly  easy  to  measure  oq  vs  tem¬ 
perature.  A  wealth  of  data  on  electrical  conductivity  has  been  amassed  for  most  metals 
and  alloys.  Thus  the  free-electron  model  is  currently  enjoying  a  great  deal  of  attention  as 
a  way  of  providing  reflectivity-vs-temperature  information  in  the  study  of  laser  effects. 

There  is,  of  course,  one  problem  in  using  cjq(T)  data  to  predict  R,  and  that  is  the 
parameter  m*/N.  It  turns  out  that  R  is  fairly  insensitive  to  this  parameter  at  infrared 
wavelengths.  To  see  this  we  show  here  some  numerical  illustration.  Define 

0  =  m*lm0 

where  mo  is  the  free  eleciron  mass,  and  then  the  parameter  /3/N  is  equivalent  to  m*/N. 
Figure  5  shows  a  plot  of  li/N)  as  a  function  of  (i/N  for  X  =  10  pm  and  various  values  of 
oq.  Here  oq  is  in  unli-  ,i  reciprocal  ohm-centimeters.  Typical  values  are,  for  example, 
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<'o  ■  1():>  ii  1  cm  1  and  ,j/,V  -  10  -  !  cm:l  for  aluminum.  Then  the  value  of  0/f/if((J/ N) 
is  ;.hout  7.2X  It)20  t'Pfl  If  vvi1  lake  a  l(Y/<  error  in  jj/N  wo  got 

A(fl//V)  -  7.2  X  102‘J  X  10' 24 

<>W/\ ) 

AR  =  0.00072. 

Since  for  those  values  R  -  0. 97300,  the  change  in  R  is  only  about.  0.1%.  We  can  obtain 
quae  good  predictions  by  the  Drudo-Lorc ntz  model  using  the  experimental  values  of  Oq 
and  the  most  simple  choice  for  gi/tV,  namely  one  free  electron  for  each  valence  electron 
per  atom  in  the  metal,  and  jj  =  1.  For  alloys,  it  is  sufficient  tc  choose  the  major  con¬ 
stituent  of  the  alloy.  For  example,  with  stainless  stc-el  we  choose  iron,  or  two  electrons 
pei  atom,  to  compute  N  and  hence  ji/.V. 


ji/H  (Cin-'l 


KiH-  3  -Sensitivity  of  R  lo  th<  parameter 
P'.V  (Ref  :i) 

Figures  (>  and  7  shew  the  predictions  of  the  free-f  iectron  theory  for  a  variety  of 
metals  and  some  compari-'on  to  ecperimental  data  |3I.  '1'he  abrupt  change  when  the 
metal  melts  is  caused  by  the  abrupt  change  in  the  d  conductivity.  Notice  in  the  com¬ 
parison  to  data  that  aluminum  films  give  values  clo.-es*  to  the  theory.  This  is  probably 
because  they  present  the  best  surfaces.  Defects,  oxide  layers,  etc.,  tend  to  trap  the  inci¬ 
dent  radiation  and  cause  the  real  surface  to  absorb  mote  radiation  than  the  ideal  surface. 
These  graphs  are  m  terms  of  absorptanee,  which  is  the  experimentally  measured  quantity, 
and.  since  metals  au;  opaque,  tl  ■  =  1  -  ‘K  ,  which  is  correct  for  specular  reflection  at  nor¬ 
mal  incidence  from  an  opaque  substance. 
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Ktg.  6  -Temporal un'  (l<'per'd>'.i<  <  of  l in'  ahsorpiimcc  al  1  0.6  V 
for  aluminum  a  ut  stainless  stool 


Lot.  us  return  to  the  expressions  for  n  and  li  to  look  at  some  limiting  forms  and  thus 
show  how  these  complete  expressions  reduce  to  simple  relationships.  Remember  that  11 
( Eq .  (2.3)  l;  determined  by  n  and  h  (Kqs.  (27),  which  are  in  turn  obtained  from  the  do 
conductivity  and  m*/<V  ( Kqs.  (28).  The  variation  of  n  and  it  with  wavelength  is  shown  in 
Kig.  8  for  a  typical  good  conductor  like  aluminum  or  copper  at  room  temperature.  Note 
that  at  long  wavelengths  n  -  k.  We  can  derive  this  by  using  Kqs.  (2b)  for  «j  and  o 2  and 
noting  that  <-.s  w  -<■  o,  ti]  -»  and  00  -*  Oqu>t.  By  substituting  these  into  Kqs.  (27)  for 
n  and  k,  we  can  readily  snow  that 


n  --  k 


r~W 


V  2  c  p  to 


(30) 


This  is  called  the  Hagen- Ruben-  lim.t.  Note  that  n  is  very  large.  Under  these  conditions 
algebra  can  be  used  to  reduce  Eq.  (23;  to 
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1 


n  +  1 


or 


1  - 

u; 


R 


2 

ti 
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Fig.  7  — Free  electron  theory  predictions  of  absorptivity  of  several  metals  at  10.6  /i. 
The  open  symbols  indicate  the  molten  state. 


WAVE  LENGTH  (M'CRONS) 


Fig.  8— n  and  h  as  functions  of 
wavelength  (Ref.  3) 
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and  Kq.  (.‘10)  can  be  .substituted  for  n  to  ni  l 


R  =  1 


2  \! 


/2to(  0 


<’o 


(31) 


This  is  tin  Ham'll- Rubens  reflectivity. 

We  can  also  comment  on  the  skin  depth.  We  have,  at  long  wavelengths  (w'o), 

X  /  2<  tjco 
O  — —  .  i  — —  • 

2.7  v  «0 


This  can  be  rewritten  as 


h 


7 


2 _ 

p0(i()  vj 


(32) 


Kquation  (32)  is  the  common  expressit  n  foi  skin  depth  used  at  long  wavelengths. 

Finally,  we  see  from  Fig.  9  that  n  and  k  reconverge  at  short  wavelengths.  I  his  is 
called  the  plasma  resonance.  To  see  this,  oi  •  must  look  at  the  behavior  of  n  and  over  a 
larger  spectrum.  We  have  already  discussed  the  long-wavelength  limiting  behavior  of  n  and 
ft.  This  is  the  Hagen-Uubens  region,  where  n  -  It.  At  short  wavelengths,  it  is  easy  to  show 
from  Eqs.  (28 1  that 


N-e  1 

0  1  t  >  •, 

in  ■  0( jio- 

AV2 
m*  co 

WAVE  LENGTH  (  MiCRONS ) 


Fig.  and  k  as  fund  mils  of  wavdon^th 
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Thus  Eq.  (27)  can  be  written,  for  large  cj,  as 


n  2 


1  - 


Ne2 

6q  m*cu2 


k2  =  0. 


(33a) 

(33b) 


Now  the  plasma  frequency  is  usually  defined  from  Eq.  (33a)  by  setting  n  =  0  to  yield 


and  thus 


Ne2 

e0m* 


(34a) 


(34b) 


We  see,  then,  that  at  very  high  frequencies  the  free-electron  model  predicts  a  transparent 
behavior  (k  =  0)  and  the  index  of  refraction  approaches  that  of  a  vacuum.  The  transitio 
to  this  transparent  behavior  takes  place  at  the  plasma  frequency,  and  it  is  a  fairly  abrupt 
transition,  as  Fig.  9  shows.  In  fact,  some  texts  call  this  transition  the  ultraviolet  catas¬ 
trophe.”  Note  that  at  to  near  Eqs.  (34)  and  (33)  are  not  valid.  For  these  frequencies 
we  must  use  the  full  expression.  If  we  use  again  the  values  of  oo  =  105  to-1  cm-1  and 
P/N  =  10"23  cm3,  which  are  appropriate  to  a  good  conductor  like  aluminum  at  room 
temperature,  the  reflectivity  looks  like  Fig.  10.  One  can  see  that,  in  terms  of  the  reflec¬ 
tivity,  the  transition  is  very  abrupt,  indeed. 


X  ifi) 


Fig.  10—  R  as  a  function  of  wavelength 
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The  optical  reflectivity  of  real  metals  is,  as  we  have  seen,  in  reasonable  accord  with 
the  free-eleciron  model  at  wavelengths  in  the  infrared.  The  surface,  however,  must  be 
nearly  perfect  for  the  predicted  reflectivities  to  be  achieved,  and,  of  course,  as  the  wave¬ 
lengths  approach  the  visible  region  band  effects  become  important  and  the  reflectivity 
shows  rapid  :! actuations  with  frequency.  The  absorptance  of  a  practical  metal  surface  is 
still  largely  a,i  empirical  matter.  For  high-power,  continuous-wave  radiation  by  a  CO2 
laser,  some  data  are  available,  but  very  little  information  on  absorptivity  as  a  function  of 
surface  temperature  under  these  conditions  is  available.  Shown  in  Table  1  are  room- 
temperature  absorptances  for  a  few  materials. 


Table  1 

Room-Temperature  Absorbtances  of  Aerospace 
Metals  and  Alloys  at  10.6  pm  for  Various  Surface 
Conditions  and  at  Normal  Incidence 


Metal  or 

Surface  Condition 

Alloy 

Ideal 

Polished 

As-Received 

Sandblasted 

A1 

0.03  3 

0.030 

0.04 

±0.02 

0.115 

±0.015 

Au 

1 

0.01 

0.02 

0.14 

Cu 

1 

0.016 

0.06 

Ag 

0.005 

0.011 

2024  A1 

0.033 

0.07 

±0.02 

0.25 

304  Stainless  steei 

0.11 

0.4 

±0.2 

Ti  Alloy 
(6A1,  4V) 

0.65 

±0.2 

Mg  Alloy 
Az-31B 

0.06 

±0.03 

Data  on  the  reflectivity  of  a  metal  during  actual  irradiation  by  a  laser  beam  is  quite 
difficult  to  obtain,  although  this  information  is  central  to  the  problem  of  laser- material 
interaction.  One  classic  experiment  along  these  lines  was  carried  out  by  Bonch-Bruevich, 
Imas,  Romanov,  Libenson,  and  Mal’tsev  in  Russia  in  1967  [4].  They  surrounded  their 
specimens  with  a  sphere  to  monitor  the  reflected  radiation,  as  shown  schematically  in 
Fig.  11.  The  output  of  the  photodetector  is  proportional  to  the  reflectance  of  the  spec¬ 
imen.  Some  of  their  results  for  steel  and  copper  are  shown  in  Fig.  12.  The  laser  pulse 
(Nd:  glass  laser,  1.06  pm),  with  a  peak  power  density  of  the  order  of  108W/cm2,  is 
shown  as  a  broken  line.  As  time  passes,  of  course,  the  laser  pulse  heats  the  surface  and 
the  reflectance  decreases.  An  especially  interesting  feature  of  these  data  is  the  shoulder. 
The  author  has  suggested  that  this  leveling  off  is  associated  with  the  surface  reaching  the 
melting  point  and  pausing  at  that  temperature  while  the  thickness  of  the  molten  layer 
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Fig.  1 2— Reflectance  of  steel  ar.u  copper  during 
irradiation  by  a  laser  beam  (from  Ref.  5;  copy¬ 
right  1969,  Clarendon  Press,  Oxford,  England. 
Used  by  permission  of  J.C.  Jaeger). 
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propagates  into  the  solid.  In  short  order,  however,  the  molten  layer  begins  to  heat  up 
and  the  reflectance  continues  to  decrease.  As  the  power  density  of  the  laser  pulse  reaches 
its  peak  and  begins  to  fall,  the  surface  temperature  can  no  longer  be  maintained,  and  as 
the  surface  cools  the  reflectance  begins  to  increase  again. 


3.  THERMAL  RESPONSE 


3.1.  Introduction 


One  of  the  most  important  effects  of  intense  laser  irradiation  is  the  conversion  of 
the  optical  energy  in  the  beam  into  thermal  energy  in  the  material.  This  is  the  basis  of 
many  applications  of  lasers,  such  as  welding  and  cutting.  We  shall  summarize  here  this 
thermal  response.  It  is  basically  a  classical  problem,  namely  heat  flow.  In  the  usual  man¬ 
ner,  we  shall  seek  solutions  to  the  equation  which  governs  the  flow  of  heat,  namely 


peg.®  if 

3 1  cx 


3T 

3x 


3  T 
3y  3y 


+  ~  Ik 


+  -r  I* 

GZ 


3  T 
dz 


+  A  . 


(35) 


We  use  here  p  for  the  density,  C  for  the  specific  heat,  T  for  temperature,  t  for  time,  and 
K  for  thermal  conductivity.  A  is  the  heat  produced  per  unit  volume  per  unit  of  time. 

In  Eq.  (35),  p,  C,  and  K  are  considered  functions  of  both  position  and  temperature,  and 
A  is  a  function  of  both  position  and  time.  In  effect,  the  equation  is  a  simple  statement 
that  the  rate  at  which  heat  accumulates  in  an  elemental  volume  clxdyclz  is  equal  to  the 
net  flow  of  heat  across  the  faces  of  that  volume  plus  the  rate  at  which  heat  is  produced 
within  the  volume. 


Thus  thermal  response  studies  consist  essentially  of  two  parts.  First,  one  needs  to 
know  the  rate  and  source  of  production  of  heat  by  the  laser,  which  yields  A.  Then  one 
solves  Eq.  (35)  subject  to  the  boundary  conditions  of  the  situation  of  interest.  This  can 
be  a  very  elaborate  task  and  frequently  can  be  done  only  with  the  aid  of  a  computer. 
There  is  a  great  deal  of  effort  among  workers  in  the  field  of  laser  effects  to  develop  an 
all-inclusive  computer  program  to  solve  Eq.  (35)  for  every  possible  situation.  However, 
the  solution  to  Eq.  (35)  can  be  no  better  than  the  knowledge  of  A,  and,  as  we  shall  see 
in  later  sections,  it  is  often  very  difficult  to  establish  A  with  any  precision  in  a  laser- 
material  interaction  situation. 


3.2.  No  Phase  Change— Semi-Infinite  Solid 


Let  us  consider  first  the  most  simple  situation.  Let  the  laser  beam  be  perfectly  uni¬ 
form  over  an  extremely  large  area,  so  that  we  have  a  one-dimensional  situation.  Assume 
also  that  the  material  parameters  are  temperature-independent  and  that  the  solid  is  uni¬ 
form  and  isotropic  and  of  semi-infinite  extent  (Fig.  13).  Finally,  assume  that  there  is  no 
phase  change;  the  rate  at  which  energy  enters  the  material  is  not  sufficient  to  induce 
melting  or  vaporization. 


First  rewrite  Eq.  (35),  using  the  fact  that  p,  C,  and  K  are  constant: 

32T  _  1_  3T  =  A_ 
dz2  *3 1  K' 


(36) 
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UNIFORM  LASER 
BEAM 


Fig.  13— Uniform  irradiation  of  a  semi-infinite  solid 


Here  we  have  introduced  x  =  K/(pC),  which  is  the  thermal  diffusivity.  Let  us  adopt  the 
convention  that  T  is  measured  with  respect  to  the  initial  (or  ambient)  temperature  of  the 
material.  This  is  possible  because  Eqs.  (35)  and  (36)  define  T  only  to  within  an  additive 
constant.  Then  we  have  as  a  boundary  condition  that  T  -*■  0  as  z  -*■  The  boundary 
condition  on  the  front  face  ( z  =  0)  depends  on  what  we  assume  for  radiative  and  con¬ 
vective  losses.  It  can  be  shown  that,  for  most  cases  of  interest,  the  rate  at  which  the 
laser  creates  heat  at  the  interface  is  overwhelmingly  larger  than  convective  and  radiation 
losses,  so  we  ignore  them  for  the  present  calculation.  Thus  the  boundary  condition  is 
that  there  is  no  heat  flux  at  z  =  0,  that  is. 


K 


ar\ 

H-o 


=  o. 


Now  consider  A.  Denote  by  I  the  power  density  of  the  laser  radiation  at  the  sur¬ 
face;  the  dimensions  of  I  are  power  per  unit  area.  The  power  density  of  the  radiation 
transmitted  to  the  surface  is  /(I  -  jR).  Then  the  power  density  as  a  function  of  z  is 

F  =  (1  -  .‘j>)/e'4jr*2/X  .  (37) 

This  follows  from  the  fact  that  the  energy  in  the  electromagnetic  wave  goes  as  £2.  Now 
to  get  the  power  transferred  per  unit  volume,  consider  elemental  volumes  of  length  dz 
and  unit  area: 


A  =  -  =  (1  -  ft)/  e-4rrfe2/X  _  (38) 

OZ  A 

The  minus  sign  appears  because  dF/dz  is  the  power  per  unit  volume  lost  by  the  radiation 
and  A  is  the  power  per  unit  volume  absorbed  by  the  material.  Finally,  we  define  the 
absorption  coefficient 


a  = 


4nk 
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which  is,  of  course,  1/5,  the  skin  depth.  Thus 


A(z,  t )  =  (1  -  5()I(t)ae-Q!:  (40) 

where  we  have  included  the  possibility  of  /  varying  with  time. 

So  the  equation  to  be  solved  is 

32T  1  3T  (1  -R)I(t)ae-az 

- K -  ,41) 

In  keeping  with  our  assumption  of  temperature-independent  thermal  parameters,  we  as¬ 
sume  further  that  R  is  independent  of  temperature.  Equation  (40)  is  valid  for  temperature 
dependent  !R  and  can  be  used  to  give  A(z,  t ,  T). 

For  metals,  or  is  a  fairly  large  number.  As  we  saw  in  Sec.  2,  k  is  of  the  order  of  100 
at  X  =  10  pm,  so  that  a  is  of  the  order  of  106  cm-1.  Hence  the  absorption  occurs  in  a 
very  narrow  layer  at  the  surface.  It  then  becomes  more  convenient  to  seek  solutions  of 


a2r  i  3 r 

az2  k  at 


(42) 


subject  to  the  boundary  conditions  that  T  =  0  at  z  -  but  with  a  specified  flux  into  the 
surface  at  z  =  0,  i.e.. 


or,  with  the  definition 


-K 


ar\ 

a*/4=o 


(!-$)/«) 


F(i)  =  (1-  «)/(«), 

~K  If)  =  F(f)> 

dzJt-  0 


(43) 

(44) 


First  examine  the  case  of  F  =  Fq,  a  constant.  This  is  appropriate  to  irradiation  by  a  con¬ 
tinuous  laser,  given  temperature-independent  material  properties.  We  note  here  only  the 
solution,  for  many  excellent  texts  on  heat  conduction  can  be  consulted  for  the  details  [5], 

The  solution  to  this  problem  is 

2 Fq  \/~Kt  _ 

T(a,  f)  = - zz -  lerfc  [a/(2  v/xt)]  (45) 


or 


T(z,t)  =  --  {y/r  e~*2l4Kt  ~i  erfc  [z/(2  %/kOj)  . 
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The  functions  which  appear  here  are  error  functions,  and  it  is  useful  to  summarize  some 
of  their  properties  and  definitions.  (See  Ref.  5,  Ch.  II.) 

The  error  function  is 


2  <*  x  .,2 

erf  ( x )  -  — ™  |  e  L  dfc 

V7T  J0 


erf  (o)  -  0,  erf  (°«)  =  1,  erf  (-a:)  -erf  (a). 


The  complementary  error  function  is 


erfc  (a)  =  1  -  erf(Aj  -  — —  /  e~'?  dY. 

v  rr  1 


The  integral  of  t lie  complementary  error  function  is 


ierfc(.T)  =  J  erfc  ((.')  dk 


or 

1  2 

ierfc(A)  - — —  e'x  -  a  erfc  (a-) 

V  T. 


ierfc(-v)  -  -|rt’  X'1  -  x  +  a  erf  (a). 


Some  derivatives  are  useful: 

b  erf  (A  |  i!  erfc  (A)  2 

()a  Oa  s/n 

<J2  erf  (A)  d-  erfc  (a)  /\ 

Pa2  Oa2  v'tt 


i)  lerfr  (a  | 

- : - -  -e.-fc(.v) 

dv 


<)2ierfc(A)  2 

- - -  — —  e  ' 

fjA2  \JU 


Now  we  can  show  that  the  boundary  condition  is  satisfied.  Using  the  first  form  of  Kip 
(4.0)  yields 


b T 
bz 


erfc  \zfd2  v^nj 
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and  siii.'C  erfe  (o)  -  1, 


jvr)  _  *0 

-0*1-0  ~  ~'K' 

□iu>  can  also  show  'hat  Kq.  (45)  satisfies  Kq.  (42). 

U'e  can  us  Kq.  (45)  to  show  what  the  front  surface  temperature  behavior  is,  under 
constant  irradiation,  by  setting  z  ~  0  so  that 


T(o,  t) 


(4b) 


fis  an  illustration,  let  us  calculate  the  time  required  to  raise  aluminum  to  its  melting 
point  for  a  power  density  of  5  kW/cm- : 

K  -  2.3  \V/cm- 

k  5=  0.9  cm2  Is 

Tm(l|l  -  600 C 

T  =  Tm(, |,  -  7’r()„m  *  600°C 

Fa  ~  (1  -  )  / 

(1  -  a',)  -  0.04 

Fu  ■=  0.04  X  5X  10:t  -  200  W/cm-. 

Then  t  -  (irK~T2  )/( 4  F,/k  )  yields  <  =  42  s.  In  practice,  it  is  very  difficult  to  melt  extremely 
thick  slabs  of  aluminum  with  even  a  high-power  laser,  as  these  calculations  suggest. 

Equation  (45).  although  derived  for  a  very  simple  case,  describes  man;,  very  important 
features  of  thermal  response  to  lasers.  First  we  shall  define  the  diffusion  length,  which  is 
useful  in  that  it  permits  a  wide  variety  of  order-of-magni!.ude  calculations  to  be  made.  The 
thermal  diffusion  length  I)  is  defined  a.^ 

D  2  y/icf  .  i  i7 ) 


Strictly  speaking,  the  thermal  diffusion  length  is  defined  as  the  distance  required  for  the 
temperature  to  diup  to  \  jc  of  its  initial  value  and  depends  somewhat  on  the  geometry  and 
the  boundary  condition.  For  most  purposes  it  is  sufficient  simply  to  take  it  as  defined 
by  Kq.  (17).  booking  at  our  solution  (Kq.  (  15)).  for  example,  we  see  that 


l>.  1 1 


2I-  o  (7<j  } 

K  V  TJ  (• 


n 
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erfe  ( 1  1 1 


26 


NRL  REPORT  77  28 


T(l),  l) 


2*’«i  |  /k7  1 
A'  l  V  7?  e 


D_ 

2 


( 0 . 1  fj  7  3 ) 


1 


Prom  Kq.  (4o), 


Thus 


T(l),  t)  =  Tin,  t) 


0.1573  s/rr 


TOT  l )  '  0.00  T(o,  7)  (48) 

in  this  case,  whereas  (1/e)  7(0.  f)  ^  0.37  7’(o,  ().  Referring  to  our  example  of  irradiating 
aluminum  for  42  s  to  reach  the  inciting  point,  we  note  that  the  diffusion  length  at  that 
time  is  given  hy 


I)  =  2  v/0.9  /  42  -  12cm, 

and,  hy  hq.  (48),  the  temperature  at  distance  O  into  the  material  is  0.09/  600,  or  about 
50cC  above  ambient. 

Now  let  us  illustrate  the  solution  (Kq.  (45))  graphically  (tig.  14).  For  convenience 
rewrite  the  equation  by  introducing  l)  -  2  \/k7  and  by  reducing  it  to  the  error  function 
erf,  so  that 


T{Z,  t  i 


2  T0  [  j) 
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Fij*.  11  —Thr  *»rrnr  function 
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(■19a) 


(19b) 


Finally,  we  define  a  dimensionless  temperature  0  -  TK/(h(jD),  so  that 

U  -  terfc  ( 7/ )  .  (49c) 

Thus  the  plot  of  the  integral  of  the  complementary  error  function  here  is  the  graph  of 
the  solution  to  the  problem  of  constant  heat  flux  on  the  surface  of  a  semi-infinite  solid. 


Now,  although  the  graph  of  Eq.  (49c)  (see  Fig.  15)  represents  very  succinctly  the 
solution  to  our  problem,  it  does  not  really  snow  how  the  temperature  varies  as  a  function 
of  position  and  time.  For  this  purpose  it  is  useful  to  look  at  the  temperature  profiles  for 
various  times  and  see  how  the  profile  changes  with  tune.  These  curves  can  be  generated 
quickly  from  II  -  ierfc  (p)  by  recalling  the  de.itiilions  of  U  and  j;  and  writing  them  in  the 
foliowing  form. 


(50b) 
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Thus,  at  a  given  time  the  d  -  ierfc  (tj)  curve  scales  according  to  Eqs.  (50);  the  basic  shape 
of  the  curve  is  unchanged,  but  it  is  stretched  one  way  or  the  other  depending  on  the 
parameters,  and  this  stretching  progresses  in  time  as  \/T.  The  case  of  aluminum  is  shown 
in  Fig.  16. 


Fig.  1 6—  Laser-induced  temperature  rise  in  aluminum  as  a  function  of  depth 


We  can  also  look  at  the  variation  of  temperature  with  time  at  a  fixed  position.  The 
variation  ol  the  surface  (z  =  0)  is  simply  T(o,  t)  ^  \/77  as  Eq.  (46)  shows.  Wherever 
z/(2  y/TTt)  is  very  small,  the  temperature  variation  will  approach  \/T  Thus,  at  any  posi¬ 
tion  T  =  y/r at  sufficiently  large  t.  The  temperature-vs-time  profiles  at  fixed  position  for 
times  such  that  z/(2>/k£)  is  not  small  can  be  calculated,  of  course,  from  Eq.  (45).  Some 
results  for  aluminum,  with  the  parameters  used  above,  are  shown  in  Fig.  17.  Notice  that 
at  z  =  10  cm  the  temperature  profile  is  far  from  the  “long-time,”  or  \/T,  behavior  even  at 
40  or  50  s,  whereas  the  surface  has  already  begun  to  melt. 

We  now  turn  to  some  order-of-magnitude  arguments.  One  such  argument  can  be 
used  to  estimate  the  power— pulse  length  combination  which  might  be  expected  to  yield 
surface  vaporization.  Consider  a  laser  pulse  that  has  the  simple  time  behavior  shown  in 
Fig.  18  and  uniformly  irradiates  the  surface  of  the  material.  The  pulse  length  is  tp  and 
the  intensity  is  such  that,  combined  with  the  reflectance,  the  absorbed  power  density  is 
Fq.  Again  we  assume  that  the  optical  energy  is  absorbed  in  a  very  thin  layer  at  the  sur¬ 
face.  Let  Dp  be  the  diffusion  length  associated  with  the  time  tp.  The  question  is 
whether  a  significant  amount  of  surface  vaporization  will  occur  before  the  pulse  ends. 

One  approach  would  be  to  use  Eq.  (46)  to  calculate  the  surface  temperature  at  the  time 
tp  and  compare  it  to  this  vaporization  temperature.  However,  this  would  ignore  the  in¬ 
fluence  of  the  latent  heats  of  melting  and  of  vaporization,  which  have  an  important  in¬ 
fluence.  We  shall  discuss  thermal  flow  with  phase  changes  later.  For  the  present  purpose 
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FiR.  1  7  —  Laser-induced  temperature  rise  on  front  surface  of 
aluminum  as  a  function  of  time 


spatially  uniform  laser  pulse, 

POWER  DENSITY  ABSORBED  F0 
PULSE  LENGTH  =  1p 


Fig.  18— Irradiation  of  a  semi-infinite  solid  by  a  pulse 
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we  can  include  them  by  considering  the  energy  required  to  melt  and  vaporize  a  portion 
of  the  material.  The  key  is  to  estimate  what  thickness  of  the  material  is  involved,  and  in 
this  order-of-magnitude  argument  we  simply  use  the  thermal  diffusion  length  for  this 
thickness.  Thus  we  set  the  criterion  for  vaporization  as 

Fptp 

r\  ^  p[Dg(Tm  ~  To)  +  Lm  +  C$>(7b  -  Tm )  +  Lv\ 

Up 


where  p  is  density  of  the  material,  Cs  and  C?  are  the  specific  heats  of  the  solid  and 
liquid,  respectively,  Tm  is  the  melting  point,  Tf,  is  the  boiling  point,  and  Lm  and  Lv  are 
the  heats  of  melting  and  vaporization.  Notice  we  are  ignoring  differences  between  the 
solid  and  liquid  for  density  and  conductivity,  as  is  appropriate  in  this  crude  argument.  If 
numerical  values  are  checked,  Lv  dominates  the  expression  on  the  right  side  of  the  inequal¬ 
ity.  For  example,  for  aluminum  Lv  =  10,875  J/g,  whereas  all  the  other  terms  contribute 
a  total  of  3,046  J/g.  Since  the  argument  is  crude,  then,  one  usually  takes 


Fptp 

Dp 


pLv 


as  the  criterion  for  vaporization  by  a  pulse.  Since  Dp  -  2  V Ktp\  we  have 


Fq  >  approx 


2  \[k  Lvp 


Some  calculations  based  on  Eq.  (51)  are  shown  in  Fig.  19.  Most  metals  fall  in  the  band 
indicated.  For  a  given  pulse  time,  at  power  densities  greater  than  the  band  indicates, 
vaporization  effects  would  be  expected  to  be  important.  Some  useful  thermal  constants 
are  included  in  Table  2. 


PULSE 

THE  BAND  INDICATES  WHERE  Eq.  (51)  LIES  FOR  MOST  METALS 

Fig.  19  — Power  density —pulse  time  criterion  for  vaporization 
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In  deriving  Eq.  (51),  we  have  been  seeking  the  power  density  required,  at  a  given 
pulse  length,  for  a  thermal  layer  to  be  vaporized.  The  same  expression,  of  course,  tells 
us  the  pulse  time  at  which  vaporization  becomes  important  for  a  fixed  power  density. 
Rewriting  Eq.  (51)  gh's  for  this  time 


K2Tlvir 
tD  >  approx  - 5 — 

4  F0~k 


(52) 


Let  us  compare  this  to  the  time  required  for  surface  vaporization  to  begin.  We  do  this 
by  using  our  solution  for  heat  flow  in  the  semi-infinite  solid  for  the  surface  (Eq.  (46)) 
and  solving  for  the  time  at  which  the  front  surface  reaches  the  vaporization  temperature: 


or 


T. 


vap 


2  Fq  ap" 
IT  V  n 


1  vap 


4F02k 


Thus,  at  tp  =  fvap  this  calculation  would  predict  that  vaporization  at  the  surface  begins. 
For  example,  at  F0  =  106  W/cm2,  vaporization  begins  at  tp  10*5  to  10"6  s,  depend¬ 
ing  on  the  metal.  On  the  other  hand,  for  a  thermal  layer  to  be  evaporated  requires,  ac¬ 
cording  to  Eq.  (52),  tp  =  10'3  s.  It  turns  out  that  both  estimates  are  useful.  In  a  later 
section  we  shall  discuss  features  of  a  more  correct  treatment,  which  accounts  for  both 
the  heat  of  melting  and  the  heat  of  vaporization  in  the  dynamic  situation  of  propagating 
solid— liquid  and  liquid— vapor  interfaces. 


3.3.  No  Phase  Change— Slab  of  finite  Thickness 

Let  us  turn  now  to  a  treatment  of  another  geometry  which  can  be  useful  in  practical 
cases,  namely  irradiation  of  one  surface  of  a  sheet  or  slab  of  finite  thickness.  Let  the  slab 
be  taken  as  infinite  in  extent  in  the  x  and  y  direction,  and  let  the  laser  irradiation  be  uni¬ 
form  over  the  entire  surface  z  =  0.  Thus  we  again  have  a  one-dimensional  situation,  as 
shown  in  Fig.  20.  The  thickness  of  the  sheet  is  taken  as  £,  the  absorbed  power  density 
as  a  function  of  time  is  F(t),  and  we  again  assume  that  the  radiation  is  absorbed  in  a  very 
narrow  layer  at  the  front  surface.  The  equation  we  wish  to  solve  is,  then, 

9 IT  _  l_  dT  = 
dz2  K  9 t 


-*  81. 


=  F(t) 


K  21 

K  9  z 


0. 


z=C 
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with  the  boundary  conditions 
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The  second  boundary  condition  states  that  -.he  rear  surface  is  insulated.  We  shall  look  at 
the  consequences  of  this  assumption  a  little  later. 

As  we  showed  for  the  semi-infinite  slab,  the  solutions  turn  out  to  be  elaborate. 
Turning  to  the  special  case  of  T(l)  =  To,  a  constant,  the  solution  is 

n7r(C  -  2)]  | 

~~r~  J  r 

'  (53) 

We  can  check  that  this  satisfies  the  boundary  conditions: 


Similarly,  the  thermal  diffusion  equation  is  satisfied,  as  the  reader  can  verify. 

Let  us  look  briefly  at  this  solution.  It  consists  of  a  linear  term  in  t ,  together  with  a 
“correcting  term,”  which  can  be  plotted  as  shown  in  Fig.  21.  In  other  words,  what  is 
plotted  is  the  term 
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U  0  2  0  4  0  6  0.8  1.0 


NOTE  that  r<  i<  l.  THE  NUMBERS  ON  THE  CURVES  ARE  VALUES  OF  Kf  ). 1 1 
F.|’.  21  —  i  ns  ,'i  fund  ion  of  1  z!\  (Ref.  f>,  p.  113) 


RfV-?)2  -  C'2  o  (  If' 

11  "  6V2"_'~  V  L*  ~ 

/» -1 


•;  t  ,iu'Z  '  n"  IF-  ■  <  I  I 
r  KnWti't*  cOS  - I. 


L  V  J 


Let  us  examine  .•  >me  special  cases.  Al  z  -  0,  for  example, 

T(>  n  ,M(  +  MTi  .  _2  V 

j  U>.  ^  <  4  Lr  I  ■>  ..2  Z—  „2  ‘  I 


A'  3  if-  i— .  »- 

/i-l 


We  can  rewrite  tins  as 


7  (<'-  1  '  "av"  7  +  A  '* ' 

Now  f  at  a  fixed  value  of  £  is  a  function  of  ll  ']  r:  Kl /V.'-,  we  can  write 

/’0V 

7V\  O  --  (?)  +  J  2  =  o<»/))- 


(54) 


(55) 
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Figure  22  shows  how  j'!  depends  on  rj,  for  2  =  0,  and  was  taken  from  the  previous  graph 
of  S'  vs  (1  -  z/E)  (Fig,  21).  Note  that  at  small  77,  i.e.,  at  Kt  «  C2,  lT  ~  0,  so  that  the 
front  surface  initially  heats  up  linearly  with  time,  as 


T(o,  f) 


£oi i  , 

KV  • 


(56) 


Fig.  Tl—  f  at  z  -  0  as  a  function  of  17 


For  large  t  values,  or  Kt  »  C2,  j\  =  o  approaches  a  limiting  value  of  about  0.33.  Thus  at 
long  times 


T(o,  t)  =  +  0.33^  ,  Kt  »  ?2.  (57) 

Here  again  we  see  linear  behavior,  but  this  time  there  is  an  additive  constant.  If  we  have 
a  very  thick  slab,  we  should  get  the  same  result  as  our  previous  solution  for  the  time  to 
reach  600° C  on  the  front  surface  of  aluminum  with  an  absorbed  power  density  of  200 
VV/cm2.  It  turns  out  that  the  limiting  form  of  Eq.  (56)  is  not  correct  because  it  ignores  the 
behavior  of  'j\  =  o(b)  at  small  77.  It  is  necessary  to  use  the  full  expression.  Thus 

KT 

v  +  -iwn)  =  ^- 

Assume  that  C  =  100  cm,  since  we  know  from  our  infinite-slab  solution  that  the  diffusion 
distance  is  12  cm  at  T  =  600° C  on  the  front.  Thus 

t)  Jz=0  (m  200X  100  °-0b9- 

Readmg  very  roughly  from  the  graph  of  D2  =  q  vs  77  gives 


36 


NRL  REPORT  7728 


r?  ft  0.004  at  ®  ft  0.065. 


Thus  our  solution  is 


»  »  °-004  =  ff  • 

which  gives,  since  K  ft  0.9  cm2/s,  a  time  of  about  44  s,  in  reasonable  agreement  with  the 
semi-infinite-slab  solution. 


Now  let  us  turn  to  a  consideration  of  the  rear  surface  temperature.  For  this  case, 
z  =  C,  so  Eq.  (53)  becomes 


T(fi,  0 


Ki 


Fq£ 

K 


o°  /_  1 

i  _  JL  N  1  i _ L.  p-Kn2n2t/<l2 

6  7T2  —  n2 

n=l 


or,  introducing  r)  as  before,  we  have 


T(2,t) 


F0i 

K 


(v  +  r?)) . 


Comments  could  be  made  here  for  the  rear  surface  temperature,  and  they  would  be  simi¬ 
lar  to  those  we  made  for  the  front  surface  temperature.  It  is  interesting  to  compare  the 
front  surface  temperature  to  the  back  surface  temperature.  This  has  a  simple  form  for 
thin  sheets,  where  Kt/22  »  1.  By  referring  to  the  graph  (Fig.  21)  of  ID  vs  (1  -  z/2),  one 
can  read  off  values  for  J)z  =  c(°°)  and  iD^oI00),  and  thus 

F0i 

T(o ,  f)  -  T(2,  t)  0.5  ~~~ 
for 


Kt/fi2  »  1. 

Notice  in  Fig.  21  that  the  limiting  values  are  approached  rapidly;  they  are  nearly  realized 
by  the  time  Kt/9.2  =1.  As  a  numerical  illustration,  if  we  have  0.3-cm-thick  aluminum, 

T(o,  t)  -  T(fi,  0  =  13°C 

For  the  same  numbers  we  used  above.  This  situation,  with  the  two  surfaces  heating  at 
the  same  rate  but  separated  by  13°C,  would  start  at  a  time  of  the  order  of  t  ft  22//c  ft 
0.1  s.  At  this  time  the  front  surface  temperature  is  about  35°C. 

Let  us  turn  now  to  a  different  sort  of  heat  input.  So  far  we  have  been  discussing 
continuous  irradiation.  Another  simple  case,  which  is  a  reasonable  approximation  under 
certain  conditions,  is  that  of  a  laser  pulse  which  is  short  enough  to  be  treated  as  a  delta 
function.  Take  again  a  slab  of  thickness  fi,  and  assume  that  the  energy  is  deposited  in  a 
very  thin  layer  near  the  surface.  F  refers,  as  before,  to  the  fraction  absorbed  by  the  ma¬ 
terial.  The  laser  power  density,  that  is,  must  be  multiplied  by  the  optical  absorptance.  In 
this  case  we  solve  the  thermal  diffusion  equation  subject  to  the  boundary  condition  that 
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-K 


a t 

d  z 


=  -if 
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'2  =  0 


=  0, 


with  the  stipulation  that  there  is  an  instantaneous  release  of  E q  units  of  energy  per  unit 
area  in  the  plane  z  -  0  at  time  zero.  This  type  of  problem  is  discussed  in  Carslaw  and 
Jaeger  [5]  and  is  most  easily  solved  by  Laplace  transform  methods.  For  our  present 
purpose  we  quote  the  solution 


*  Q* 

‘■S' 

C*}  *>s 

II 

•♦o 

sii 

[l  *  2  f;  cos(nlz)  c-«»2»2</e2l 

n -1  '  ' 

(58) 


In  Eq.  (58)  we  have  introduced  Eq,  the  energy  per  unit  area  in  the  pulse.  Thus, 


E  o 


F(t)  dt . 


For  the  case  under  consideration,  F(t)  is  considered  to  be  a  delta  function. 

Equation  (58)  is  the  basis  for  a  scheme  used  quite  frequently  for  the  measurement  of 
thermal  parameters  [6],  This  scheme  consists  of  using  a  thin  sheet  of  the  material  to  be 
studied  and  irradiating  uniformly  one  surface  with  a  very  short  laser  pulse  while  monitor¬ 
ing  the  temperature  rise  induced  on  the  back  surface.  If  one  knows  Eq  ,  and  if  the  as¬ 
sumptions  of  no  heat  loss  are  valid,  the  experiment  can  yield  values  of  both  specific  heat 
and  thermal  conductivity.  One  adjusts  the  pulse  energy,  and  hence  Eq,  so  that  the  in¬ 
duced  temperature  rise  is  small.  In  this  way  the  values  of  specific  heat  and  thermal  con¬ 
ductivity  are  representative  of  essentially  the  ambient  temperature  of  the  material. 

To  see  how  this  is  applied,  rewrite  Eq.  (58)  for  the  back  surface,  z  =  (2: 


T(Z,  t ) 


Eqk 

KZ 


1  +  2  ^  (~l)n  e"  <n2*2tlZ2 

n*l 


(59) 


If  we  introduce  a  characteristic  time  tc  =  C2/X7r2,  Eq.  (59)  looks  approximately  like  the 
curve  shown  in  Fig.  23.  Here  we  have  also  introduced  a  characteristic  temperature  Tc  ~ 
Eqk/KZ  and  plotted  T/Tc  vs  t/tc,  or 


X 

Tc 


1  +  2  )  (-l)ne'n2«/'t  . 

L — j 

n=  1 


(60) 


Essentially  the  experiment  consists  of  monitoring  the  temperature  as  a  function  of  time 
and  fitting  it  to  Eq.  (60).  This  can  be  done  quite  readily.  First,  the  long-term  tempera¬ 
ture  rise  Tx  yields  the  specific  heat  because 


T 
1  c 


1 
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,/fC 


Fij*.  23  —  Normalized  back  surface  t^mpnralur#*  response  to 
a  delta  function  heat  pulse 


and,  on  substituting  for  Tc, 


or,  since  k  ■  Ki/jC , 


K»k  . 

KH  ' 


<ai ) 


7 hi.-.  technique  of  measuring  specific  heat  is,  of  course,  not  unique  to  pulsed  lasers  It  is 
sometimes  referred  to  as  the  slat/  calorimeter.  The  accuracy  of  the  method  depends  on 
knowing  A'o,  which  is  frequently  difficult  to  ascertain  with  laser  radiation.  In  same  ap¬ 
plications  Eq.  ihf)  is  used  to  calculate  /•.(-,  the  energy  actually  absorbed  from  the  laser 
pulse,  by  using  materials  o  known  specific  heat. 


The  pulsed  laser  measurement  technique  is  especially  suited  to  determining  thermal 
diffusivity.  The  magnitude  of  the  back  surface  temperature  rise  depends,  as  we  saw,  on 
the  energy  which  is  coupled  into  the  material,  and  the  may  be  difficult  to  know  with 
any  precision.  However  the  time  dependence  of  the  hack  surface  temperature  is  inde 
pendent  of  the  energy  input  and  is  controlled  only  by  the  diffusivity  k.  A  simple  way 
to  derive  k  from  a  temperature-time  profile  can  be  seen  From  Fig.  23.  One  measures  the 
time  required  for  the  measured  temperature  response  to  reach  some  fraction,  say  one- 
half,  of  its  limiting  value.  Let  us  call  this  time  t\p>.  It  can  be  shown  numerically  [6| 
from  Kq.  (60)  that 


rjn 
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J.  T.  SCHRIEMPK 


fl/2/7c  -  137. 


Using  tc  ■■  k-jiKn'1)  yields 


.  j_37 
n  2  ll  n 


(62) 


Thus  a  measurement  of  together  with  the  thickness  of  the  specimen,  immediately 

gives  the  thermal  diffusivity.  If  one  knows  An  this  experiment  gives  values  of  both  the 
specific  heat  and  the  diffusivity.  and  hence  if  one  knows  the  density  p  the  experiment 
gives  the  thermal  conductivity. 

This  technique  lias  been  applied  often  at  rather  high  temperatures,  usually  in  the 
lOOCUC  range  and  above.  At  these  temperatures  steady-state  methods  of  measuring  the 
thermal  conductivity  are  difficult  to  apply  because  radiation  losses  are  so  large.  In  the 
laser  flash  technique  the  radiation  loss  goes  like  7'^  Tn' ,  whore  To  is  the  starting  or 
ambient  temperature,  established  by,  say,  a  furnace,  and  is  defined  above.  This 
radiation  loss  can  he  made  quite  small  by  adjusting  A’o  so  that  is  only  a  few  degrees 
larger  than  'i  j  Since  the  precise  value  of  A‘()  is  difficult  to  establish,  these  experiments 
typically  measure  only  the  thermal  diffusivity,  not  the  conductivity 

A  final  remark  on  the  criterion  for  the  applicability  of  Eq.  (60)  to  slab  heatn  g  con¬ 
cerns  the  limits  on  the  laser  pulse  duration.  No  laser  pulse  is,  of  course,  a  true  b  function. 
Our  solution  would  be  expected  to  be  correct  for  laser  pulse  times  which  are  short  com¬ 
pared  to  the  time  it  takes  the  hack  surface  to  respond.  The  response  times  are  of  the 
order  of  t, ,  so  we  have  the  criterion 


or 

Ip  «  t'-'/(K7T7). 

Calculations  that  include  explicitly  the  time  dependence  of  the  laser  pulse  |7]  indicate 
that  our  6-fu  .ction  solution  in  error  by  less  than  about  27,  provider’  that  (,,  is  less  than 
ot  equal  to  about  4  V.  of  tc.  Some  typical  value;:  of  tc  with  specimens  1  mm  thick  are 
given  below. 


Aluminum 
Stainless  steel 


k  tc  -  V'-'/KTr2 

(em2/s)  (ms) 


0.85  1.2 

0.0523  19 


Typical  laser  pulse  length  ~  are  about  1/2  to  1  ms  for  the  so-called  “normal  mode” 
lasers,  and  thus  with  1-mm-thiek  rtecimens  tin*  technique  would  be  fairly  accurate  for 
stainless  steel  hut  not  very  good  for  aluminum.  Thicker  specimens  would  help,  hut  tins 
would  make  the  rear  surface  temperature  rise  smaller.  If  our  laser  pulse  has,  say,  20 
.I/cm2  in  it  and  we  use  the  1 0.6  pm  absorptance  quoted  earlier  for  as-received  surfaces, 
the  anticipated  temperature  rises  at  the  hack  of  the  j-inni  specimens  would  he  as  shown 
beiow. 
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A'o 

P 

C 

?oo  -  kuipCSi 

(1 

(J/crn-) 

(g/cm2) 

(cC) 

Aluminum 

0.0  1 

0.8 

2.7 

1.05 

2.8 

Stainless  steel 

0.4 

8.0 

8.0 

0.628 

16 

We  see  Oie  need  to  use  shorter  laser  pulses  hut  with  the  same  amount  of  energy.  Another 
solution  would  he  to  carry  out  a  more  detailed  heat-flow  calculation.  Both  tailoring  of 
the  pulse  shape  and  more  detailed  calculations  are  usually  employed  in  current  applica¬ 
tions  of  laser  flash  techniques  |8|. 


3.4.  Melting 

Consider  first  the  case  of  a  semi-infinite  slab  melted,  with  instantaneous  melt  re¬ 
moval,  as  indicated  in  Fig.  24. 


F 


o 


Z  =  0  1=  Ut 


SEMI -INFINITE 
MEDIUM 


Fill.  2  1  --  Irradiat  inn  of  a  semi  i n rin i t r  slab  with 
in.star.l uneous  inell  removal 


We  solve  this  by  considering  ourselves  moving  along  with  the  interface.  But  to  do 
this  we  have  io  reconsider  our  heat-flow  equation,  which  is 


h-'V  J_  '22’ 
x  i)t 


0  . 


Recall  that  this  was  derived  by  noting  th  -  rale  ai  which  heat  accumulated  in  an  elemental 
volume: 


0  . 


If  the  medium  is  moving,  an.  additional  amoun:  of  heat  pCT  is  flowing  in  at  rate  V. 
So 
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..  37’  .  ,,  bT  , 

~h  - —  becomes  -a  -t—  +  VpCT. 

3a  az 


Thus,  for  moving  media,  the  heat-flow  equation  is 


pC  —•  ■*  ~~  f  K  —  '  +  VpCT]  -=  0, 
dt  3a  \  dz  j 

provided  we  are  not  generating  any  heat  in  the  solid. 

Thus 


37’  327’  ,,  37’ 

—  K - —  +  V  - 

dt  3  z2  6z 


0. 


(63) 


So  we  can  find  at  least  a  steady-state  solution  to  fie  problem  by  basing  our  coordinate 
system  at  the  interface  and  letting  the  material  flow  in  at  some  rate  to  be  determined. 

(.'all  the  rate  -{t,  where  U  is  a  positive  number.  The  solution  will  be 

7’  =  Tm  e  <■’*''*  (641 

where  tc  is  diffusivity  of  the  solid,  Tm  is  the  melting  temperature,  and  z‘  is  the  distance 
from  the  melting  front.  That  is,  this  is  the  solution  to 


3T  <)-T 


=■-  0, 


as  we  can  verify.  First,  bT/bt  -  0  because  this  is  tire  steady-state  profile,  or 


K 


bz"2 


3_ 
3  a 


0. 


(65) 


Of  course  T  -  T,„  at  z‘  -  0,  and  T  -*  0  as  z'  -*  !t  is  trivial  to  show,  by  differentia¬ 

tion,  that  Kq.  (64)  is  the  solution  to  Kq.  (65). 

do  calculate  U,  we  use  energy  balance.  I\)  must  raise  the  material  to  Tm ,  its  melting 
point  (Tm°C  atiove  ambient),  and  then  melt  it.  'thus  in  tune  ill  the  energy  put  into 
thickness  .Aa'  (where  U  ~  *1 z'/iXt )  must  be  given  by 


or 


F^J 

<Aa 


Lp  +  cr,„p 


!■„  -  p\I-  +  CTm\  i; . 


(66) 


Thus  we  ran  write,  for  the  steady-state  temperature  profile, 
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T 


>0 _  I 

/'{L  +  | 


(67) 


where  Fq!{/j(L  +  C7’m)|  is  the  velocity  of  the  melting  front.  This  solution  also  would  be 
appropriate  to  sublimation,  where  L  is  then  the  heat  of  sublimation  and  Tm  the  sublima¬ 
tion  temperature.  Note  that  this  is  the  semi-infinite-slab  approximation  and  cannot  be 
used  to  estimate  the  time  to  penetrate  a  slab  of  given  thickness. 

Let  us  consider  aluminum,  with  Fq  ■■■  200  VV/cnvL  Taking  a  more  accurate  value  of 
the  molting  point  than  in  earlier  examples,  T,„  ■  f. 40° C  (the  actual  melting  point  of 
fi60J(l  minus  room  temperature  of  20°C).  If  we  put  in  the  other  values  of  the  param¬ 
eters,  T  falls  off  as  shown  in  Fig.  25. 


j'  (Cm) 


Kit;.  -Temperature  profile  in  aluminum  wilh 
instantaneous  melt  removal 


As  another  illustration,  consider  plexiglass.  Plexiglass  is  rapidly  eroded  hy  10.6-pm 
laser  radiation,  by  a  process  that  is  essentially  sublimation.  Since  it  couples  extremely 
well  (cl  '  1.0)  and  has  a  very  low  thermal  diffusivity  (x  =?  ]  0  1  ein2/s),  it  can  he  used  to 
make  “burn  patterns”  of  the  beam.  That  is,  the  depth  of  the  erosion  at  a  given  point  is 
linearly  proportional  to  the  energy  density  incident  at  that  point.  We  cun  understand 
this  by  applying  heat-flow  concepts.  Let  us  apply  Kq.  (67)  and  interpret  L  as  the  heat  of 
erosion.  Since  C  1.1  J/g-C  and  l'm  *  200  "C  and  L  *  1000  J/g,  for  a  crude  ca  cula- 
tion  we  can  ignore  CTm.  The  density  is  about  1.1  g/em:i.  We  then  have 

T  -  200  c  ',/oJ  . 

Since  (1^1,  for  a  typical  [rower  density  like  o  kW/cm*  we  have  f'0  =  5  X  1(),!,  so  that 

T  ~  200  e  ‘  r,x  O’4-2' 

Hence  the  temperature  profile  is  confined  to  an  extremely  narrow  region  near  the  eroding 
surface.  The  rate  of  erosion  is,  hy  Kq.  (Ob)  with  (,T.:I  ignored. 
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(/’  =>  4.5  cm/s  . 

Now  we  can  see  why  plexiglass  is  useful  for  monitoring  the  beam  profile,  and  why 
it.  can  pick  up  fairly  fine  structure  in  the  beam.  After  irradiation  for  a  time  t  the  pene¬ 
tration  depth  is  Ui  -  (F’0n/(p/A  and  this  should  be  large  compared  to  a  thermal  diffusion 
length  l)  -  2  \^Kl  if  the  pattern  is  to  r<«.-«al  fine  structure.  Otherwise  thermal  diffusion 
would  “wash  out’’  the  pattern  by  distributing  the  energy  in  a  radial  direction.  Thus 

_ _  f*Q  t 

2  %/kT  <<  — —  ■ 

pL 

For  the  numbers  we  used  above  at  t  =  1/4  s,  we  have  a  depth  of  1  cm.  Thus 


;  i 

2  cl  O' 4  X  —  ==  0.1  mm  <<  1  cm 
V  4 

and  we  see  the  criterion  is  well  satisfied 

Let  us  look  at  a  more  complete  problem,  namely  melting  by  laser  radiation  of  a 
slab  of  material.  One  basic  problem  is  what  happens  to  the  melted  material.  (We  will 
ignore  the  vaporization  question  for  now.)  There  are  two  cases  which  are  fairly  amenable 
to  numerical  solution.  They  are  the  “fully  retained  liquid"  case,  in  which  all  the  liquid 
is  presumed  io  stay  in  place,  and  the  “full  ablation"  case,  in  which  the  melt  is  presumed 
to  disappear  magically  as  soon  as  it  forms.  The  latter  ease  might  correspond  to  the  pres¬ 
ence  of  a  heavy  windstream  which  blows  away  the  melt. 

Lo<  king  first  at  the  full  ablation  e;  e,  we  have  the  situation  shown  in  Fig.  2b.  7-j 
is  the  temperature  in  the  solid,  and  the  iront  surface  at  z  -  0  first  warms  up  to  the  melt¬ 
ing  point  'm  (above  ambient)  and  then  begins  to  move  to  the  right.  .S  denotes  its  posi¬ 
tion  as  a  function  of  time.  When  .S’  -  U,  the  process  iS  over,  and  we  call  this  time  tf.  We 
denote  by  tm  the  time  at  which  the  front  surface  begins  to  melt.  The  field  equation  is 


z  -0 


T  =  Tm 


\ 


7 


INSULATED  BACK  SURFACE 


z;S  j  l--£-  z 
SOLID  REGION,  T  =  T?  (  z,t  ! 
Fi*>.  2i'y  Fully  ablated  cast* 
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The  above  boundary  conditions  are  nonlinear,  and  a  solution  in  analytical  form  is 
very  difficult.  This  is  due  to  the  presence  of  the  moving  boundary  and  appears  in  the 
second  boundary  condition,  which  states  that  the  boundary  moves  at  a  rate  dS/dl  deter¬ 
mined  by  a  balance  between  the  heat  of  melting  L,  the  heat  input  Fq,  and  the  heat  flow 
by  thermal  conduction. 

One  relationship  must  hold  for  this  problem;  it  follows  from  energy  balance.  The 
total  energy  put  in  per  unit  area  is  Fotf ,  and,  since  the  material  is  simply  heated  to  Tm 
and  melted,  this  energy  goes  solely  to  those  processes.  Thus 

Fotf  =  p(L  +  CTm  ) . 

This  is  convenient,  for  one  can  check  numerical  solutions.  More  important,  it  gives  a 
first-order  estimate  of  the  time  needed  to  melt  through  materials  by  laser  radiation. 

Turning  now  to  the  fully  retained  liquid  case,  we  have  the  following  set  of  equations. 
The  definitions  are  the  same  as  above,  except  that  subscript  1  now  refers  to  the  molten 
state,  whereas  subscript  2  is  still  the  solid  state  (see  Fig.  27).  The  field  equations  are 

a2r2  i  ar2 

U  '  «1  IT  =  0  o«(<v  (solid) 

82r,  1  3Ti 

IT-  "  ~i  IT  =  0  '*><•<  ‘r  (liquid) . 

The  boundary  conditions  are 
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T  =  Tm 


INSULATED  BACK 
SURFACE 


1  =  0 


i=  S 


z  =  J. 


LIQUID  REGION, T  =  T,  (  z,t)  \ 

SOLID  REGION,  T=T2  (r,t) 


FiR.  27  — Fully  retained  liquid  case 


-  K2 

0T2\ 

**  to 

=  ^0 

0  <  t 

-  Kl 

<>zLo 

=  *b 

tm  < 

t  <  tf 

*  W-,,, 

-  *2 

dT<>\ 

=  ~PL 

dS 

dt 

t,n  * 

/A 

k2 

ar2\ 

=  0 

0  <  f 

<  lf 

^2)z  =  5  ^1^  =  5  '^n 


<  t 


tf- 


The  initial  conditions  are 


T2(z,  o)  -  T20 


0. 


Rather  than  discuss  this  problem  in  detail,  we  pass  on  to  the  more  practical,  although 
more  complex,  case  of  vaporization.  For  the  fully  retained  liquid  case,  suffice  it  to  say 
that  the  retained  liquid  has  a  shielding  effect,  and  thin  causes  the  time  to  reach  melting  at 
the  back  surface  to  be  longer  than  in  the  ablated  model.  Some  typical  values  of  melt- 
through  time  for  0.2-cm-thick  material  with  Fq  ~  2  k\V/cm2  are  given  below. 
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Ablated 

Retained 

(s) 

(s) 

Aluminum 

0.32 

0.37 

Stainless  steel 

•1.0 

4.5 

3.5.  Melting  and  Vaporization 

We  present  here  without  derivation  some  results  for  the  ease  of  a  slab  of  material, 
insulated  on  the  surfaces,  subjected  to  undorm  and  continuous  irradiation  [9J.  These 
are  one-dimensional  calculations.  It  is  assumed  that  the  melt  is  fully  retained  until  it 
reaches  the  vaporization  temperature,  where  it  disappears.  Then  we  have  the  case  illus¬ 
trated  in  Fig.  2H. 


LIQUID  I  SOLID 


INSULATED  BACK 
SURFACE 


Fit;.  2H  Mi'loni!  and  vn|>oii/.,ition,  with  fully  ■‘••turned  liquid 


Here  .5'2  is  the  position  of  the  liquid-vapor  interface  and  Si  the  position  of  the  solid- 
liquid  interface.  This  problem  his  been  solved  numerically  at  NRL  19).  and  we  shall 
show  some  results.  The  assumpt’ons  are  that  in  each  phase  the  thermal  properties  are 
independent  of  temperature,  in  the  curves,  the  following  definitions  are  used: 

F0V 

o Q  =  — ; - 

b-Ks„i„t 

Ksolid  f 

7  '-nr- 

V  (1  i  m  *  so  I  ill  it* 
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L  -  latent  heat  of  fusion 

subscript  0  ■  ambient  temperature 

so  =  LoiL ,  with  Lv  the  latent  heat  of  vaporization. 

In  these  equations  T  is  understood  to  he  in  degrees  Celsius  and  represents  Lire  actual  tem¬ 
perature.  Although  the  thermal  conductivities  of  liquid  and  solid  are  allowed  to  differ, 
the  specific  heats  are  assumed  to  be  the  same. 

Note  that  when  Fq  -*  0  and  Fo  -*•  ■*>,  we  get  certain  easy  limits.  Fo*-  Fq  (1  no  va¬ 
porization  can  take  place,  the  melting  is  small,  and  we  approach  tne  fully  ablated  limit. 
On  the  other  hand,  as  Fb  ~y  °°,  all  the  liquid  should  he  vaporized  by  (/•,  the  time  the  back 
surface  melts. 

So  in  this  limit 

~  F?i  F  +  (  T,n  ~  7o)Cs(,|,d  +  (  I  v  ~  T  m )  ^"liquid  +  ]  ( t  ij  ,J°)  , 

whereas 

*otf  -  P? \L  +  {Tm  ~  '/b)C'so|„i]  (F0  ~  0). 

The  limiting  values  are  re.  >nted  hy  t h<»  asyn’ptotes  of  the  curves  (Figs.  29-41) 
indicated  by  dashed  lines  on  the  plots  of  ixQ  vs  t.  Note  on  these  plots  that  the  dasred 
lines  ar<-  at  45“ ,  or  have  a  slope  of  -1.  These  are  log-log  plots,  so  the  asymptotes  can  be 
described  by 

log  («Q)  -  log  0  -  log  Tf 

where  c’  -*  c' x  as  c/Q  -■»  <»,  and  c'  c  <j  as  aQ  -*  zero.  If  we  take  the  antilog, 

o  O  -  c jy  , 

substituting  in  the  definition  of  a Q  and  t f  gives 

Fo?  ,  l! 2 

—t — —  =  t  - — 

/'F^solid  Ksolnlf/ 

or 

Fotf  •  P?U  ’  . 

By  comparison  with  the  Fq  -*  0  and  the  Fu  °°  limits  above  we  can  see  that 
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Fig.  31— Power  density  vs  burn-through  time  for 
aluminum  alloy  (Ref.  9) 


The  numerical  values  of  the  thermal  parameters  which  were  used  in  generating  these  solu 
tions  are  included  as  a  separate  table  (Table  3)  in  addition  to  the  graphs  (Figs.  29-34). 


Fig.  32  — Rear  surface  temperature  rise  for  stainless  steel  (Ref.  9) 
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Provided  one  knows  Fq,  these  solutions  are  reasonable  estimates  for  the  time  to 
penetrate  a  metal  specimen  with  a  laser  beam.  In  application,  however,  one  must  con¬ 
sider  the  actual  size  of  the  beam.  These  solutions  will  be  useful  for  effects  in  the  center 
of  the  beam  if  the  diffusion  length  is  small  compared  to  the  beam  radius,  or  if  for  times 
up  to  and  including  the  melt-through  time  tf  the  beam  radius  R  is 

R  >  Isfictf. 

In  terms  of  the  parameter  Tf  this  becomes,  upon  squaring  both  sides, 


4.  EFFECTS  OF  PULSED  LASER  RADIATION 
4.1.  Power  Levels  of  Pulsed  Lasers 

Highly  intense  pulses  of  short  duration  can  be  produced  in  a  variety  of  ways.  Typi¬ 
cally  it  is  done  by  creating  a  large  population  inversion  by  the  injection  of  electrical 
energy  from  the  discharge  of  large,  highly  charged  capacitors.  In  these  systems  all  the 
energy  is  produced  in  a  burst,  the  duration  of  which  can  be  made  quite  short.  We  shall 
not  discuss  the  various  techniques  by  which  these  pulses  are  created,  but  in  Table  4  we 
simply  note  some  commonly  obtained  values  [2]. 


Table  4 


Laser  Type 

Pulse 

Length 

Power 

Energy 
Per  Pulse 

Ruby  (normal  mode) 
Ruby  (Q-switched) 
Ruby  (mode-locked) 
C02  TEA 

C02  e-beam 

C02  shock  tube  gdl 

0.1-1  ms 
10'8  s 

10’11  s 

10X  10'6  s 
2  X  10"5 

3 X  10‘4  s 

10-100  kW 
1-10  GW 
0.1-1  TW 

100  MW 

50  MW 

0.3  MW 

1-50  J 

1-1C  J 
0.1-1  J 

100  J 

1000  J 

100  J 

From  the  above  table  it  is  apparent  that  with  beam  areas  of  the  order  of  1  cm2  extremely 
high  power  densities  can  be  obtained,  and,  although  the  pulse  lengths  are  short,  the  total 
energy  in  each  pulse  is  considerable.  The  available  power  densities  range  as  high  as  1012 
W/cm2. 

Practically  speaking,  one  is  usually  interested  only  in  power  densities  below  the 
breakdown  threshold  of  air  because  at  higher  power  densities  the  energy  never  reaches 
the  target.  These  breakdown  levels  are  functions  of  wavelength,  spot  size,  and  pulse 
length,  and  depend  as  well  on  the  contaminants  in  the  air.  Typical  values  are  109  W/cm2 
in  “clean”  air  at  STP  for  C02  laser  pulses  with  duration  of  about  10'6  s  and  longer.  At 
shorter  pulse  lengths  the  threshold  is  somewhat  higher,  becoming  1010  W/cm2  at  10‘8  s 
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and  1011  VV/cm2  at  10" 10  s.  In  the  infrared  region,  the  breakdown  threshold  scales  with 
the  square  of  the  frequency. 


4.2.  Material  Vaporization  Effects 

We  shall  first  discuss  the  effect  of  high-power-density  laser  pulses  on  materials  from 
the  point  of  view  of  target  vaporization,  and  shall  assume  that  the  vaporizing  surface  is 
not  shielded  from  the  radiation  by  the  vapor.  In  this  case  we  can  show  that,  in  addition 
to  thermal  input  to  the  target,  there  is  a  strong  pressure  built  up  on  the  target  surface 
due  to  recoil  from  the  blowoff  of  the  vapor.  The  integral  of  this  pressure  over  the  time 
of  the  laser  pulse  imparts  a  net  impulse  to  the  target.  There  arises  then  the  possibility  of 
inducing  stresses  large  enough  to  create  gross  mechanical  changes,  such  as  spall  and  de¬ 
formation,  by  pulsed  laser  irradiation. 

To  calculate  the  pressure  applied  to  a  surface  by  a  laser  pulse,  we  start  with  a  con¬ 
sideration  of  the  vaporization  process.  We  use  a  one-dimensional  calculation  because  in 
most  cases  of  interest  the  beam  radius  R  is  larger  than  the  thermal  diffusion  length  during 
the  pulse  time  tp ,  or 


R  >>  2  \/  k  tp  . 

We  shall  avoid  consideration  of  thin  targets,  so  that  C  is  also  large  compared  to  the  diffu¬ 
sion  length.  In  this  case  we  can  calculate  the  time  tf,  required  for  the  front  surface  to 
reach  the  vaporization  temperature  Tv  from  the  semi-infinite  slab  result  of  Eq.  (46), 
which  is 


or 


or,  since 


'r 
4  u 


2Fq  fKtb 

K  V  7T 


tb 


7 tK2T* 

4  FqK 


K 


K 

PC' 


(b 


KpCtf 


(68) 


In  applying  Eq.  (46)  in  this  way,  we  ignore  the  molten  layer  and  assume  that  the  values 
of  K ,  p,  and  c  appropriate  to  the  solid  can  be  used.  This  is  not  as  gross  an  approxima¬ 
tion  as  it  may  seem,  because  at  these  power  densities  the  molten  layer  is  very  thin. 

Once  the  material  on  the  surface  reaches  the  boiling  point,  the  surface  begins  to  erode 
at  a  rate  t/R  given  by  energy  consideration,  as  we  saw  in  section  3. 
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Ps [C'solid  +  t‘m  +  ^liquid(^i'  T'm)  +  ^v] 

To  simplify  the  calculation  we  take  6'^^  =  Csojjd  =  C  and  ignore  Lm  by  comparison  to 
Lv.  Then 


Us  = 


Fq 

Ps(Uv  +  CTV) 


(69) 


Here  we  have  used  ps  for  the  density  of  the  solid.  So  after  the  time  tj,  given  by  Eq.  (68) 
the  surface  begins  to  evaporate,  and  it  recedes  at  the  rate  Us.  By  conservation  of  momen¬ 
tum  it  must  be  true  that 


PvUv  =  psUs  (70) 

where  pv  and  Uv  designate  the  density  and  velocity,  respectively,  of  the  evaporation  prod¬ 
ucts.  Thus  we  have 


PvUv 


Fq 

Lv  +  CTV 


(71) 


by  combining  Eqs.  (69)  and  (70). 

To  see  how  the  pressure  exerted  on  the  surface  is  related  to  density  and  velocity, 
note  that  the  pressure  on  the  surface  is  just  the  pressure  of  the  evaporation  products.  To 
calculate  this  pressure,  consider  particles  which  move  a  distance  A z  in  time  A(  under  the 
pressure  P  and  thereby  acquire  a  velocity  V.  The  pressure  (force  per  unit  area)  must 
equal  the  rate  of  charge  of  momentum  (per  unit  area)  so  that 

(pAz)  V 
At  ' 

That  is,  pAz  is  the  mass  per  unit  area  which  is  brought  to  velocity  V  in  time  At  by  the 
force  per  unit  area  P.  V’  =  Az/At,  and  so  P  =  pV2.  Thus,  in  our  specific  case  of  density 
pv  and  velocity  Uv  there  is  an  associated  pressure,  given  by 

F  =  PvUt2  .  (72) 

We  could  compute  the  pressure  from  this  expression  if  we  knew  pv  and  Uv.  How¬ 
ever,  we  only  know  the  product  pvUv,  from  Eq.  (71).  We  need  another  relationship, 
which  we  simply  take  from  the  ideal  gas  law, 


where  R  is  the  gas  constant  and  A  the  molecular  weight.  Denote  R/A  by  Cf '  and  use 
Eq.  (72),  then 
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or 


Ut.uf  -  PvK'Tu 

IK  =  777;. 


C’3) 


Upon  combining  Kqs.  (71)  through  (73),  we  get  the  desired  relation!, hip. 


*b  v^X 

i.,i  +  67; 


Since  the  specific  heat  of  metal  is  typically  Hi i A  we  can  approximate  c" '  by  (1/3)  C  to 
yield 


_1 _ -T>  7  CTV 

7.5  |iu  +  CT„]  ' 


(74) 


Finally,  we  compute  the  sjx’eific  impulse  delivered  during  the  pulse,  wii.ch  is  the  force  per 
unit  area  multiplied  by  the  time  over  which  it  acts,  and  we  gel 


or 


U>) 

,  1  Fq>/cZ  '  Ti  KpCT?\ 

‘m  '  x/3  \L„  +  C'T^I  4  f2  j 

\  u  / 


(7  5a) 


In  terms  of  the  energy  in  the  pulse,  K o  -  Fo tp,  -bis  can  be  written  as 


l  m 


1  Fu  \fCTu  /  j[  KpCl , 
x/3  \LU  *CTV]  "  4  to2  ^ 


(75b) 


A  word  about  units  is  in  order.  It  has  become  conventional  to  quote  impulse  in 
units  of  dyne-s,  and  specific  impulse  in  dyne-s/em2.  If  we  use  J/em2  for  energy  density, 
J/g-X'  for  specific  heat,  and  J/g  for  heat  of  vupomation,  we  have 


I /ml 


J/cm'2 

\/j/? 


-  cm " 2  s/J  -ti 

~  cm  " 2  7 id7  erg  g 
or 

|/,n!  "  TuT7  (dynp-s)/cm2  . 


This  unit,  dynes/cm'2,  is  called  a  tap. 


Thus,  in  taps, 
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=  lta  x  io3 


sfCT,_  /  7T  KpCTff 


\Lu  +  C7’tr] 


l1 


So" 


(76) 


Note  t.iat  Kq  (76)  predicts  a  threshold  value  of  /->  for  impulse  production  at  a  given 
pulse  length  t/;.  This  is  due  to  the  criterion  we  introduced  for  vaporization;  vaporization 
must  commence  before  the  end  of  the  pulse  or  there  will  be  no  impulse.  The  threshold 
is  give  n  by 


or 


£ 

4 


KpCTy 
c  2 


tn  = 


1 


£o)th 


Vtr 

2 


7't,  \JKpCtp  . 


This  vaporization  model  also  predicts  that,  at  very  large  /.0,  the'  impulse  per  unit  area  is 
directly  proportional  to  the  energy  density  with  a  constant  coupling  coefficient ,  given  by 


U 1 


1.83  X  1<)3 


( f-  i;  +  i'7v) 


This  is  the  iimii  at  which  vaporization  begins  essentially  instantaneously  with  respect  to 
the  pulse  length  and  vapor  products  are  produced  for  the  entire  pulse. 

Some  numerical  values  arc1  illustrated  below  and  in  Kig.  35.  f,'o  is  in  J/cm-  and  tp 
ii  /is,  so  t hat  !m  is  in  taps: 


Kc  r  titanium 


For  aluminum 


In, 


' 

8.01  K0  1 


6  23  tp/h  q 


In, 


6.94  i: o 


1  -  33.9  tp/K 


Vhc*  above  mod  •!  illustrate.;  the  principles  involved  ri  generating  impulse  by  laser 
vaporisation.  In  fact,  in  predicting  threshold  values  it  gives  results  which  are  within  a 
factor  if  two  of  experimental  measurements.  It  has  been  refined  110)  by  a  calculation 
which  ;  mounts  for  the  fact  that  7',,  is  probably  not  a  handbook  value  that  comes  from 
measure  ments  at  atmospheric  pressure,  but  rather  a  different  value'  appropriate  to  the 
dynamic  and  high-pressure  situation  created  by  the  laser  induced  vaporization.  In  this 
refinement,  T,,  is  determined  from  the  kinetic  model  of  vaporization,  which  predicts  that. 
l\  ■  c„  i  x p  |  1,,-HR'T,,)] ,  where  e0  u  the  speed  of  sound  in  the  solid  and  R'  the  gas 
constant  per  gram.  When  this  is  doin',  the  thresholds  agree  very  well  with  theory.  How¬ 
ever,  as  /-.'ii  (and  hence  /■', since  tn  is  constant)  is  increased,  experiments  show  that  de¬ 
livered  impulse  does  not  inc-rea.se  indefinitely  but  begins  to  fall  off.  This  is  due  to  the 
onset  of  absorption  of  laser  energy  by  the  vapor  produc  ts  and/or  the  heated  air  near  the 
target.  We  turn  now  to  a  consideration  of  this  problem. 
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.0  - 


Fit;.  :l5--Cou|>lin>!  coefficient  vs  energy  density  for  titanium 


4.3.  Effects  from  Absorption  of  Radiation  in  tire  Plume 

The  plume  of  vaporized  material  blown  off  the  target  becomes,  at  some  power  den¬ 
sity,  hot  enough  that  it  or  the  air  begins  to  absorb  the  laser  radiation.  The  onset  of  this 
process  is  not  thoroughly  understood,  and  the  ignition  of  these  so-called  absorption  waves 
is  the  subject  of  a  great  deal  of  research.  Proper  t-eatment  of  the  problem  depends  on, 
among  other  things,  computing  the  onset  of  ionization  and  the  rate  of  absorption  of  light 
by  the  electron:;  and  also  accounting  for  both  cascade  processes  and  relaxation  processes 
m  a  full  dynamic  sense.  We  shall  not  treat  this  problem  here.  Rather  we  shall  look  at 
some  crude  models  which  show,  in  a  semiquantitative  way,  various  features  of  the  ab¬ 
sorption  process. 

First  note  that  the  decoupling  of  the  absorption  from  the  material  surface  due  to 
shielding  by  the  plume  depends  on  the  wavelength  of  the  radiation.  Recall  from  Eqs.  (34) 
that  at  the  plasma  frequency  the  reflectivity  of  a  'Tree-electron”  metal  drops  sharply 
from  a  value  near  unity  to  essentially  zero.  If  we  assume  that  the  coupling  to  the  plume 
u  due  to  the  light  interacting  with  electrons,  Kq.  (34a)  is  valid.  Using  the  mass  of  the 
free  electron  gives 

Pp  =-  K.97  X  1U:1A,I/2  (77) 

for  the  plasma  frequency  in  hertz,  when  <V  is  in  electrons  per  cubic  centimeter.  This  can 
lie  rewritten  in  terms  of  the  corresponding  wavelength  \p  to  yield 

N  -  (1.12X  l<;::i)/v  (78) 

where  is  m  centimcti  is.  At  a  given  wavelength  the  plume  is  transparent  until  the  elec¬ 
tron  density  reaches  the  value  given  by  Kq.  (78),  where  there  will  he  a  transition  to  a 
condition  hi  which  the  plume  absorbs  and  reflects  the  radiation  and  thereby  shields  the 
material.  For  10.(i-gim  radiation,  shielding  begins  at  1  0 1  ‘-J  electrons  per  cnr!,  for 
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1  .(Hi  pm  (Nd)  at  10-'  electrons  per  nir’,  and  lor  0. ($1 '  '3  pm  (luby)  at  2.3  X  1 0  2 1  elec¬ 
trons  per  cm3. 

When  the  electron  density  reaches  a  high  enough  value,  the  beam  decouples  from 
the  surface  and  presumably  the  pressure  due  to  material  blowoff  will  drop.  To  yet  some 
idea  of  the  order  of  magnitude  of  the  energy  density  for  a  given  pulse  length  where  this 
process  begins,  let  us  simply  assume  that  cutoff  begins  when  the  front  surface  reaches 
the  temperature  at  which  the  material  is  fully  ionized.  his  should  predict  an  upper 
limit,  for  full  ionization  is  obviously  not  required.  Foi  .  sample,  solids  have  N  **  1023 
cm-3,  whereas  we  only  require,  at  10.6  pm,  N  >  1019  cm'3.  For  simplicity  assume 
that  melting  and  vaporization  processes  can  he  ignored,  and  again  use  for  the  front  sur¬ 
face  temperature  rise  the  simple  expression 

,,,  2  to  rKj 

'  “  ~r  v'v- 

A  typical  ionization  temperature  for  a  metal  would  be  about  75,000  C.  Using  simply  the 
values  of  K  and  k  for  the  solid,  we  get  for  titanium, 

F0  n/V  *  5.7  X  10  !  W  s,/2/cm2  . 

Using  K o  -  Fyt,  this  can  be  rewritten  as 

Ft-  =*  57  (79) 

where  is  understooci  to  be  time  in  microseconds. 

Figure  30  shows  some  data  taken  bv  Dr.  Rudder  of  the  Air  Force  Weapons  Labora¬ 
tory,  at  two  pulse  lengths,  1.2  and  11  p s,  with  1.06-pm  radiation  and  titanium  targets 
|11],  The  lines  marked  “Anisimov  predictions"  are  calculated  from  the  vaporization 
model  of  the  previous  section  with  the  refined  method  for  determining  7',,.  (This  was 
first  done  in  the  Soviet  Union  by  Anisimov  [19].)  The  experimental  data  agree  very  well 
at  values  of  F  near  threshold.  Note  that  Kq.  (79)  for  estimating  the  onset  of  shielding  is 
roughly  consistent  with  these  oat  a,  although  the  experimental  onset  of  shielding  is,  as  one 
might  expect,  fairly  gradual.  The  line  on  the  graph  marked  LSD  Predictions  refers  to  a 
theoretical  estimate  based  on  the  idea  that  the  laser  light,  when  it  couples  into  the  blow- 
off,  can  create  an  explosionlike  shock  wave  in  the  air  which  travels  up  the  beam,  absorb¬ 
ing  the  radiation  energy  in  the  process.  This  Iaser-supported-detonation,  or  LSD,  wave  is 
one  form  of  laser-supported-absorption  wave.  We  wdl  discuss  these  waves  next. 

Once  the  coupling  of  the  radiation  with  the  ejected  vapor  (and  perhaps  the  air) 
reaches  a  sufficient  level,  the  absorption  region  begins  to  behave  in  a  fashion  character¬ 
ized  hy  hydrodynamic  dissipation  of  the  energy  coupled  into  it.  For  now,  let  us  ignore 
the  ignition  problem.  The  absorption  region  typically  propagates  up  the  laser  beam  in  a 
way  that  is  determined  bj  the  medium  in  which  it  propagates  (usually  air)  and  also  by 
the  balance  between  the  power  being  fed  ir.  by  the  laser  and  the  relaxation  processes 
which  dissipate  the  power.  Three  types  of  laser-suppoi  ted -absorption  waves  are  usually 
identified.  'Typical  power  levels  at  which  they  appear  and  their  typical  velocities  of 
propagation  are  indicated  below  for  ll).f-pm  radiation  and  targets  in  air  at  standard  tern 
peraiurc  anti  pressure  |12). 


59 


J.  T.  SCHRIEM1' 


I'  in  36— Specific  impulse  delivt'iod  to  solid  targets  b>  I  ,06-pm 
laser  radiation  { R ef-  1 1 ) 


Type  of  Wave 

Power  Level 
of  Laser  Flux 
(VV/cm?- ) 

Velocity  of 
Propagation 
(em/s) 

Lusor-supported-detonation 

107 

105 

wave  (LSDW) 

Laser-supported  combustion 

10« 

10“ 

wave  ( 1 ,3CVY ) 

Plasmatron 

10"1 

0 

The  LSD  wave  propagates  as  a  shock  wave,  i.e.,  at  supersonic  velocity,  whereas  the 
DSC  wave  moves  more  slowly  unci  relaxes  by  thermal  conduction.  The  plasmatron  is  at 
rest,  with  the  energy  input  being  balanced  by  reradiation  and  convective  losses  into  the 
atmosphere.  Although  we  discuss  these  effects  here  in  the  section  on  pulsed  lasers,  they 
are  just  ns  valid  for  continuous  radiation.  Since  pulsed  lasers  art'  the  most  convenient 
devices  for  reaching  these  power  levels,  especially  for  LSD  waves,  absorption  waves  are 
usually  considered  under  pulsed  effects. 

Hydrodynamic  theory  can  be  applied  to  model  these  waves.  The  problem  was  first 
solved  in  the  Soviet  Union  by  Raider  |ld].  Detonation  waves  can  be  discussed  most, 
readily  because  tne  hydrodynamic  equations  reduce  (o  fairly  simple  expressions,  so  we 
shall  consider  them  in  some  detail.  A  few  remarks  about  combustion  waves  will  come 
later. 
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We  can  derive  conditions  for  the  steady-state  behavior  of  a  detonation  wave  by  con¬ 
sidering  conservation  of  mass,  momentum,  and  energy  at  the  detonation  front.  For  this 
purpose  we  do  not  concern  ourselves  with  how  the  process  starts  but  presume  that  a 
detonation  wave  has  been  formed  and  is  propagating  at  some  steady  rate  as  sketched  in 
Fig.  37.  The  absorption  region  is  propagating  to  the  right  at  a  steady  velocity  u.  We 
assume  that  it  is  very  thin  and  can  be  treated  as  a  detonation  front.  Thus  u  is  the  det¬ 
onation  velocity.  The  temperature  and  density,  etc.,  of  the  air  go  through  very  rapid 
cnanges  in  the  very  short  distance  V.  Note  that  this  wave  propagates,  in  this  treatment, 
in  air,  and  thus  our  results  will  he  independent  of  target  material. 


LASER  BEAM 
INTENSITY  I 


Fig.  37 --Temperature  and  density  profiles  typical  of  a 
laser-supported  deioiialion  wave  (Ref.  12) 


In  this  discussion  “behind  the  front”  refers  to  the  high-teniperature-and-pressure 
region  imn  »‘diately  to  the  left  of  the  absorption  region  in  Fig.  37.  “Ahead  of  the  front” 
is  to  the  right  in  the  sketch  and  refers  to  ambient  air  conditions.  Mote  we  have  given  the 
beam  a  finite  radius  li  and  thus  will  have  to  consider  lateral  expansion.  First  let  us  do 
the  oiie-dimensional  problem  and  assume  that  the  detonation  front  propagates  simply  as 
a  plane  wave. 

Behind  the  front  let  p,  P ,  and  e  be  the  density,  pressure,  and  internal  energy  per 
urit  mass,  respectively,  end  let  po.  Ai,  and  e<i  he  the  same  variables  ahead  of  the  front. 
Define  the  velocities  witl  respect  to  a  coordinate  system  moving  with  the  front  at  the 
detonation  rate  u.  Then  the  ambient  gas  •  into  the  front  with  the  speed  u,  and  we 

define  v  as  the  speed  with  which  the  high-pressure  gases  leave  the  front.  We  can  now 
write  down  t.he  conservation  equations  for  mass,  momentum,  and  energy  across  the  deto¬ 
nation  front.  These  equations  are  based  on  flow,  that  is,  they  are  in  the  terms  of  “per 
unit  area,  per  unit  time,”  The  equation  for  mass  is 


Mass 


Po“  =  PV ■ 


(HU) 
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The  conservation-of*momentum  condition  results  from  equating  the  impulse  to  the 
change  in  momentum.  Now  impulse  is  force  multiplied  by  time,  but  in  the  “per  unit 
area,  per  unit  time”  sense  this  becomes  simply  pressure.  Since  mass,  in  this  flow  con¬ 
cept,  is  pv,  momentum  is  (pii)u.  Hence  we  have 

p  -  Po  =  ~{(pv)v  -  (p0u)u)  , 

and,  if  we  ignore  Pq,  which  is  much  smaller  than  P,  we  have  for  momentum 
Momentum 


P  +  pu 2  =  pQU2 ' 


(81) 


The  conservation-of-energy  condition  follows  from  similar  considerations.  The  dif¬ 
ference  in  energy  flow  on  each  side  of  the  front  must  be  balanced  by  the  work  done  on 
the  gas  ( Pqu  -  Pv)  and  the  energy  absorbed  from  the  laser  beam,  which,  using  our  earlier 
notation,  is  F.  F  =  CT /,  where  Cf  is  the  absorptance  of  the  gas  in  the  absorption  region. 

Thus  we  have 


pu[e  +  \  v 2)  -  p0u  (e0  +  |  u2)  =  Pqu  -  Pv  +  F. 


If  we  use  P0  *  0  and  eo  0  and  substitute  from  Eq.  (80),  we  get  for  energy 


.  1  2  1  2  P  F 

e  +  V1  -  77  - - +  - 

22  p  Pqu 


or 

Energy 


e 


+ 


F 

Pqu  ' 


(82) 


Our  goal  is  to  use  these  conservation  laws  to  predict  the  pressure  P  behind  the  front 
and  ultimately  the  pressure  transmitted  to  the  target.  For  now  assume  that  F  is  known, 
and,  of  course,  the  ambient  air  density  p0  is  known.  Thus  we  have  three  equations  and 
five  unknowns,  P,  p,  v,  e,  and  u.  To  proceed  we  need  to  invoke  some  equation  of  state 
for  the  gases,  and  we  shall  simply  assume  that  the  ideal  gas  law  holds.  Thus  we  have 

P  -  pR'T  (83) 

where  R'  is  the  gas  constant  per  unit  mass,  or  R'  =  R/A,  where  A  is  the  molecular  weight. 
(Since  the  wave  is  presumed  to  be  in  air  in  this  treatment,  .4  would  be  the  average  mo¬ 
lecular  weight  of  air.  Taking  A  for  air  to  be  29.4  g/mol  gives  R'  =  2.84  X  106  erg/(g-°C) 
and  is  consistent  with  the  ideal  gas  law  and  with  p0  =  1.29 X  10"3  g/cm3  for  air  at  0°C 
and  1  atmosphere  =  106  dyne/cm2.)  Now  Eq.  (83)  essentially  introduces  another  un¬ 
known,  the  temperature  T,  so  we  need  to  add  the  expression  for  the  energy  of  an  ideal 
gas,  which  is 
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e 

where  7  is  the  ratio  of  the  specific 
take  7  =  1.4. 

Now  we  have  four  equations  (Eqs.  (80 )-( 82 ),  (84))  in  five  unknowns  (P,  p,  u,  e,  and  u). 
To  get  the  final  condition  we  use  the  criterion  for  detonation,  which  is  that  the  velocity 
of  the  high-pressure  gases  behind  the  front,  relative  to  the  front,  is  equal  to  or  greater 
than  the  local  speed  of  sound.  Intuitively  this  seems  reasonable,  for  propagation  of  shock 
waves  is,  by  definition,  in  excess  of  the  speed  of  sound.  The  criterion  can  be  properly 
derived  from  a  consideration  of  the  thermodynamics  of  the  situation,  but  we  shall  not  do 
so  here  [14].  Since  we  shall  be  interested  in  the  minimum  value  of  I  (or  F)  required  to 
sustain  a  detonation  wave,  we  take  u  equal  to  the  speed  of  sound.  For  an  ideal  gas  the 
sound  speed  is  (7P/p)1/2,  so  we  have  our  last  condition, 


R'T 


(84) 


7-1  (7  “  1)P 

heats,  7  =  Cp!Cv.  For  our  purposes  it  is  sufficient  to 


u2  = 


yP 

P 


Before  discussing  the  algebra,  let  us  collect  the  equations- 

p0u  =  pv 

X 

P  +  pu2  =  PqU2 


A  P  A  1  2  1  9  .  -F 

e  +  —  +  —  =  —  u-  +  - 

P  2  2  .  Pqu 


e  = 


v 2  = 


(7  -  l)p 

yP 


Combine  Eqs.  (80)  and  (81)  to  yield  expressions  for  u2  and  e2: 

Pp 


u2  - 


V 2  = 


P0(P  ~  P 0) 
Ppo 


P(P  -  Po) 

Now  use  Eqs.  (85)  and  (87)  to  eliminate  v2  and  p  so  that  we  can  get 


(85) 

(SO) 

(81) 

(82) 

(84) 

(85) 

(86) 

(87) 


P_  _  1  +  7 
PO  7 


(88) 


which  is  one  of  the  equations  we  need,  namely  p  in  terms  of  the  known  quantities  Po 
and  y.  Now  we  can  use  Eq.  (81)  to  get  P  in  terms  of  u  by  eliminating  u2  with  Eq.  (85)  to 
get 
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p  -  BLu2. 

1  +  7  ' 


(89) 


We  need  one  more  relation  to  complete  the  solution,  namely  u  in  terms  of  F.  This  vvi'l, 
by  Eq.  (89),  give  us  P  in  terms  of  F.  To  get  this  we  use  Eq.  (82)  and  replace  e  via  Kq. 
(84)  and  v2  via  Eq.  (85).  Thus  Eq.  (82)  becomes 

P  P  1  yP  =  _1  2+  F 

(7  -  1  )p  +  P  2  9  2  Poll  ' 

If  we  use  Eq.  (89)  to  eliminate  P,  vve  get 


or 


or 


u2  _ 1 Pi j  po 

1  +  7  [(7  1)  P  +  P 


PO  j  _1_  o 
P  \  2 


+  ¥ 
r'U11 


2  7(1  +  7), 


F 

PO" 


1  ££ 

2  P 


If  we  use  Eq.  (88)  for  p^ip,  we  get 


2 


(72  -  1) 


+ 


r 

po"  ’ 


Finally  we  arrive  at 


or 


2  F 
Po 


u 


(90) 


The  equations  which  rcprcsent  the  solution  for  the  detonation  wave,  then,  are 

_P_  =  1  +  7 
PO  7 


P 


Po"2 

1+7 


(88) 

(89) 
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U 


2(72  1  )F 

Pa 


(90) 


These  three  equations,  together  with  the  ideal  gas  law,  represent  the  formal  solution 
to  the  propagation  of  the  laser-supported-detonation  wave.  Given  the  temperature  be¬ 
hind  the  front,  and  since  Eq.  (88)  defines  p,  we  could  calculate  /'  and  hence  u  and  finally 
the  F  required  to  support  it.  However,  this  does  not  really  solve  the  problem.  What  we 
wish  to  discover  is:  given  the  laser  intensity  /,  will  an  LSD  wave  be  supported?  To  an¬ 
swer  this  question,  we  need  to  consider  the  distance  it  takes  for  the  laser  radiation  to  be 
absorbed.  We  also  need  a  more  realistic  situation  then  the  simple  plane  wave. 

First  we  note  that  the  beam  has  a  finite  radius  R  and  that  lateral  expansion  can  take 
place.  The  order  of  magnitude  of  the  radial  expansion  velocity  will  be  the  speed  of 
sound  ca.  To  maintain  the  detonation,  we  must  replace  the  energy  lost  to  expansion  by 
energy  put  into  the  absorption  region.  To  simplify,  let  us  assume  all  of  the  laser  beam 
energy  is  absorbed  in  the  distance  2.  (Actually  the  beam  intensity  only  falls  by  1/e  in 
the  distance  2.)  We  define  At  as  the  time  for  the  shock  front  to  move  a  distance  2,  or 
At  =  2/u.  In  this  time  the  radial  expansion  is  the  amount  caAt.  Now  irR2IAt  is  the 
energy  deposited  by  the  beam  in  the  cylindrical  volume  shown  in  Fig.  38  (caAt  «  R). 
But  the  energy  in  this  volume  after  expansion  is  approximately  equal  to  its  volume  multi¬ 
plied  by  its  internal  energy  poi  unit  volume.  Thus 

7T R2IAt  Poe[TrRH +'2itRcaAti)  . 


Fiy.  38  — Cylindrical  volume  which  absorbs 
beam  energy  via  expansion 


Since  At  -  2/u, 


InR2 


Poeu 


ttR2  +  2? tRc3  - 
3  u 


Note  that  for  7  =  1.4,  Eq.  (88)  predicts  that  p.  =*=  2pp  and  hence,  using  pv  =  pca  =  PqU, 
u  =  2ca.  So  we  have,  after  some  algebra, 


p0eu 


(91) 
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But  this  equation  simply  represents  the  rate  at  which  we  must  put  energy  into  the  absorp¬ 
tion  volume  in  order  to  maintain  the  conditions  we  assume  to  carry  out  our  detonation 
wave  calculation,  namely  a  plane  wave  propagating  by  absorption  of  laser  energy  in  a 
distance  (!.  Thus  the  energy  flow  per  unit  area  from  the  laser  beam  is 


F  = 


/ 

1  +  (t//?)  ' 


(92) 


Finally,  we  can  complete  the  problem  if  we  know  the  absorption  length  t.  To 
compute  C  we  need  to  invoke  some  model  of  the  ionized  air.  For  this  purpose  it  is  suf¬ 
ficient  to  assume  that  free  electrons  absorb  the  light  and  that  the  electrons  come  from 
singly  ionized  atoms.  We  shall  not  derive  the  expressions  which  we  need  but  simply 
quote  them.  There  are  two  relationships.  The  first  of  these  is  the  Saha  equation  (14a], 
which  relates  the  fraction  of  atoms  ionized  a  to  the  absolute  temperature  T  and  the  ion¬ 
ization  potential  1  of  a  single  atom: 


3/2 

O'2  O  El  21  -  ( //;'/' 

1  -  «  2  go  P  \  h  2  /  e 


(93) 


In  this  equation  m  is  the  mass  of  the  atom,  mo  is  that  of  the  electron,  k  is  Boltzmann's 
constant,  and  h  is  Planck’s  constant.  The  statistical  ^weights  of  the  ground  state  of  the 
atom  and  its  first  ionized  state  are  go  and  g\,  respectively.  Typically  g\  =  go  -  1.  In 
terms  of  known  constants,  then,  the  Saha  equation  gives  us  the  degree  of  ionization  as  a 
function  of  temperature. 

Knowing  the  degree  of  ionization,  we  can  get  the  absorption  length.  Again  we  sim¬ 
ply  quote  the  relationship  [12],  which  assumes  that  the  light  is  absorbed  by  inverse 
bremsstrahlung.  The  expression  is 


\  _4  (  2  n 

<!  3  1 3m()liT 


modlw)'*  \  m 


which,  at  the  temperatures  of  interest  and  for  10.6-qm  radiation,  becomes  (for  hv/l;T«  1 


J.  _  i.  f  \112  e6/i2  .  p2a- 

V  3  \3moj  f7i2moc(fm)-  (/r7')'5'" 

in  Kq.  (94)  v  is,  of  course,  the  frequency  of  the  laser  radiation,  c  is  the  speed  of  light,  and 
e  is  the  electronic  charge.  By  combining  Kqs.  (93)  and  (94)  we  can  calculate  V  in  terms 
of  temperature  T,  density  p,  and  known  parameters.  A  typical  value  of  i  for  air  is  of  the 
order  of  13  eV.  (For  O 2  i  is  12.1  eV,  for  No  it  is  15.f>  ('V.)  Thus  we  have 

«  -  fiP ,  T) 

or.  since  /;  1  +  y)!y  '  2 P(),  we  can  get  a  relationship  between  (’  and  T  and  hence 

between  /  and  T,  via  Kqs.  (90)  and  (92).  Typical  results  are  shown  in  Fig.  39. 
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Fig.  39—  Relationship  between  /  and  T  Tor  a  lasersupported-detonatior.  wave  (Ref. 
12).  This  plot  is  based  on  a  more  realistic  expression  for  the  equation  of  state  of 
the  gases  than  the  ideal  fias  law,  but  the  ideal  cas  law  gives  a  similar  result. 


Here  we  have  assumed  a  beam  radius  of  10  cm.  The  important  point  is  that  there  is 
a  minimum  in  the  /-vs-jT  relation.  We  identify  this  as  the  minimum  flux  lm  required  to 
maintain  an  LSD  wave.  Associated  with  it  is  the  temperature  Tm  of  the  high-pressure 
region  at  the  detonation  front.  We  can  then  use  our  detonation-wave  relationships,  Eqs. 
(SS)  through  (90),  to  get  the  pressure  behind  the  front,  or  equivalently,  get  pressure  from 
Tm  via  p  2po  and  the  ideal  gas  law. 

We  shall  turn  to  a  calculation  of  the  pressure  on  the  target  in  a  moment.  First  note 
that  the  radial  expansion  concept  imposes  a  natural  criterion  for  the  difference  between 
a  combustion  wave  and  a  detonation  wave.  The  time  for  radial  expansion  is  R/ca  whereas 
the  time  for  passage  of  the  absorption  region  is  t/u.  If  the  detonation  condition  is  to  be 
maintained,  radial  expansion  times  must  be  larger  than  propagation  times  in  order  for  the 
high-pressure  region  to  move  as  a  shock  front  and  not  dissipate  itself  radially.  Hence 
Rlca  >  V./u.  We  have  already  noted  that  u  =  2 ca  or,  crudely,  u  ~  ca,  so  that  Ii  >  fi,  or 
V.iR  <  1,  is  the  condition  for  detonation  waves.  If  t  becomes  larger  than  R  the  absorp¬ 
tion  region  is  large,  the  relaxation  in  the  radial  direction  is  important,  and  the  process 
called  a  combustion  wave  takes  place.  This  can  be  treated  in  a  similar  fashion  to  the 
detonation  wave,  but  the  hydrodynamic  equations  do  not  take  the  simple  form  of  Eqs. 

(80 1 -(82).  We  shall  not  treat  combustion  waves  in  this  report.  The  solution  in  Fig.  39 
is  for  a  detonation  wave  and  hence  is  valid  for  V/R  <  1.  The  limit  t/ff  =  1  is  shown  in 
the  figure  by  a  dashed  line. 

Finally,  compute  the  impulse  delivered  to  a  target  by  a  laser  beam  of  intensity  /  just 
sufficient  to  maintain  a  detonation  wave  [  II ).  The  beam  has  a  pulse  duration  We 
wish  to  calculate  the  effect  on  the  target  due  to  the  “explosion  products”  behind  the 
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absorption  region.  These,  of  course,  expand  in  ail  directions  and  create  a  pressure  on  the 
target.  To  demonstrate  the  effect  we  shall  use  a  very  simple  model,  namely  a  model  of 
cylindrical  expansion.  We  consider  that  the  absorption  region  has  propagated  a  distance 
Z  by  the  end  of  the  laser  pulse,  and  at  that  time  we  have  created  a  cylinder  of  high- 
pressure  gas  which  has  a  radius  equal  to  the  beam  radius  R,  a  length  Z,  and  a  pressure 
Pd  given  by  Eqs.  (89)  and  (90)  above,  with  P  =  Pd.  This  cylinder  is  then  allowed  to  ex¬ 
pand  radially  at  a  speed  estimated  to  be  the  speed  of  sound  ca.  Then  we  get  the  impulse 
delivered  to  the  target  by  integrating  the  force  on  the  target  due  to  the  pressure  in  the 
expanding  cylinder  during  the  time  the  cylinder  expands  from  R  to  the  target  radius  Rt- 
For  Rt  very  large,  the  integration  is  stopped  when  the  cylinder  pressure  drops  to  atmos¬ 
pheric  pressure.  The  model  is  sketched  in  Fig.  40  at  the  time  t  =  tp. 


TARGET 


\.Z. 


Z  =  utc 


Fig.  40— Radial  expansion  model  for  impulse  delivered  to  a  target  from  a  laser-supporled-delonation  wave 


The  model  might  be  expected  to  be  valid  if  Z  »  R  and  if  tp  «  the  time  required 
for  radial  expansion,  either  to  Rt  or  to  atmospheric  pressure,  to  take  place.  We  also  are 
assuming  that  impulse  due  to  target  vaporization  is  negligible,  i.e.,  that  the  detonation 
wave  is  formed  very  early  in  the  laser  pulse. 

We  take  the  radial  expansion  to  be  at  constant  temperature.  Then  P  times  the  vol¬ 
ume  t"  of  the  cylinder  is  a  constant.  Since  our  model  presumes  only  cylindrical  expan¬ 
sion,  we  have  the  condition  that 


Pr 2  =  constant  (95) 

where  r  is  the  radius  of  the  cylinder  and  R  <  r  <  Rt-  We  shall  need  this  relationship  in 
the  derivation  of  the  impulse.  Let  the  impulse  be  I'm,  and  let  F  be  the  force  on  the  tar¬ 
get  due  to  the  pressure.  This  gives 

rl(r=Rr) 

/;  =  f<u 

Jl(r=R) 

where  the  upper  limit  is  understood  to  be  valid  only  where  P  is  greater  than  atmospheric 
pressure  at  r  =  Rt.  Since  the  radial  expansion  rate  is  c„ 
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Since  /.  <  Rp,  we  shall  ignore  the  /?2/2  term.  The  specific  impulse  /m,  which  is  /„, 
divider  by  the  area  of  the  laser  beam  7 tR2  is  then 

,  _  m - V*  +  7 

m  2R  v*dPo  y 

Recalling  the  expression  in  Eqs.  (89)  and  (90)  for  and  ignoring  i/R  with  respect  to 
unity,  we  can  write 


h  71  +  y\  (  P o  \  2(72  -  l)  f 
IR  72  j  P°  (l  +  T )L  Po 


This  simplifies  to 


i  |2(72  -  1)) 1/3 


p?'3/1'3. 


The  expression  in  braces  is  nearly  equal  to  unity  for  typical  value  of  7  (say,  7  =  1.4),  so 


■4-pS'3'1'3. 


Now  the  energy  per  unit  area  in  the  beam  is  E 0  =  ltp ,  and  we  can  write  the  coupling 
coefficient  Im/E 0  as 


11_  2/3 //,  1  /  3  t’2/3 

IR  P°  /VP 


This  is  the  equation  of  the  straight  line  marked  LSD  Predictions  shown  in  Fig.  36,  where 
the  calculations  were  done  for  the  parameters  appropriate  to  the  1.2-ps  pulse  length. 

Several  important  consequences  of  the  LSD  wave  are  seen  in  Eq.  (97).  One  is  that 
the  coupling  coefficient  is  reduced  as  Eq  becomes  larger,  which  tells  us  that  we  cannot 
create  an  arbitrarily  large  impulse  at  a  target  by  simply  increasing  the  energy  in  the  laser 
beam.  In  fact,  when  Eq.  (97)  is  considered  to  be  correct  at  high  Eq  and  the  results  of 
the  vaporization  model  (see  Eq.  (76)  and  Fig.  36)  are  used  at  lower  values  of  Eq,  there 
is,  for  a  given  pulse  length,  an  optimum  value  of  Eq  for  transferring  the  largest  amount  of 
impulse  to  a  target.  For  the  1.2-ps  pulse  illustrated  in  Fig.  36  the  optimum  value  of  Eq 
is  about  22  J/cm2,  and  this  is  in  reasonable  accord  with  the  data.  Of  course  the  specific 
impulse  lm  per  se  goes  as  Eq^,  so  larger  Eq  will  create  larger  /,„.  However,  this  slow  in¬ 
crease  of  Im  with  Eq  is  a  very  inefficient  way  to  impart  stress  to  a  material.  A  better 
scheme,  perhaps,  would  be  to  use  multiple  pulses  at  the  optimum  Eq  value. 

Another  consequence  of  the  LSD  wave  is  a  lack  of  dependence  of  impulse  on  the 
parameters  of  the  target  material.  The  same  impulse  is  produced  independently  of  the 
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target.  This  is  in  accord  with  experiment.  When  /  is  well  into  the  lange  where  LSD’s  are 
formed,  the  measured  impulse  is  the  same  for  all  target  materials.  Some  data  taken  by 
NRL  [15]  are  shown  in  Fig.  41.  In  this  graph  we  see,  in  accord  with  the  vaporization 
model,  a  strong  dependence  of  In,  /Eq  on  material  type  at  lower  power  densities,  whereas 
at  the  high  power  densities  typical  of  LSD  formation  the  values  of  Im/E q  are  the  same  for 
all  materials.  In  this  range,  however,  a  target  area  dependence  appears.  The  target  area 
dependence  shown  is  for  aluminum.  Here  again  the  general  behavior  predicted  by  Eq. 

(97)  can  be  seen. 


AVERAGE  POWER  DENSITY  (WAT  T/CM2) 


Fig.  4  1  —Coupling  coefficient  as  a  function  of  power  density  (Ref.  1  5) 


t 


We  have  not  yet  considered  the  radius  at  which  the  expanding  cylinder  reaches  at¬ 
mospheric  pressure.  Call  this  radius  Rq.  As  explained  above  Pr2  is  constant,  so 

P0R02  =  PdR2 

where  Pq  is  atmospheric  pressure,  106  dyne/cm2.  Thus, 

P0  =  PP0'1/2Pj/2. 

From  Eqs.  (89)  and  (90),  for  Pd  this  becomes 
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,-1/2 


P  0 


2(72  -  1)/ 


n2/3 


«o » R/,o-yrT7L“ 

Again  the  factors  involving  y  are  nearly  unity,  so 

R0  =  RPompl,6lM. 
Upon  substituting  /  =  EqRp ,  we  get 

*o  =  i?Po1/2Po/64/3^1/3- 


(98) 


So  if  i?0  '-s  less  than  Rt>  one  replaces  Ry  by  R o  in  Eq.  (97).  Target  sizes  larger  than  /Jo 
will  all  receive  the  same  impulse. 

A  numerical  illustration  is  useful.  Let  Eq  =  1000  J/cm2  and  tp  =  100  ps.  Suppose 
the  beam  radius  is  1  cm.  The  power  density  I  is  about  107  W/cm2,  so  we  expect  a  LSD 
wave,  if  we  take  Ry  =  5  cm  for  the  target  radius,  Eq.  (97)  yields  (with  po  =  1.29 X  10'3 
g/cm3) 


Jm  5£ 

E0  2 


(i.29 x  io-3)2'/[(io-4)J/3ao10)2/3] 


dyne-s 

erg 


or 


'm 

Ro 


*  7X  10~7  dyne-s/erg 


—  ~  7  dyne-s/J  . 

In  this  example  expansion  to  atmospheric  pressure  would  take  place  at  a  radius  given  by 
Eq.  (98): 


R0  ---  (106r1/2(1.29X  10-3)1/6(1010)1/3(10-4r1/3 
where  we  have  used  erg/cm2  for  Eq.  Then 

R0  *  15  cm  . 

Thus  targets  with  radii  of  15  cm  and  larger  wouid  exhibit  a  maximum  coupling  coefficient 
of  ( 1 5/5 )2  times  7,  or  about  63  dyne-s/J.  If  we  compute  Im.  we  have  63,000  dyne-s/cm2, 
that  is  63,000  taps,  as  the  maximum  specific  impulse  from  this  laser  pulse.  Since  this  im¬ 
pulse  is  delivered  in  times  of  the  order  of  magnitude  of  the  laser  pulse  time,  this  corre¬ 
sponds  to  a  pressure  of  roughly  (6.3  X  1 04 >/ 1 0" 4  *  6  X  108  dyne/cm2,  or  abcut  600 
atmospheres. 
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